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ABSTRACT 

Aims. The earliest phases of massive star formation are found in cold and dense infrared dark clouds (IRDCs). Since the detection 
method of IRDCs is very sensitive to the local properties of the background emission, we present here an alternative method to search 
for high column density in the Galactic plane by using infrared extinction maps. Using this method we find clouds between 1 and 
5 kpc, of which many were missed by previous surveys. By studying the physical conditions of a subsample of these clouds, we aim 
at a better understanding of the initial conditions of massive star formation. 

Methods. We have made extinction maps of the Galactic plane based on the 3.6 - 4.5 yum color excess between the two shortest 
wavelength Spitzer IRAC bands, reaching to visual extinctions of ~ 100 mag and column densities of 9 x 10^^ cm"^. From this we 
compiled a new sample of cold and compact high extinction clouds. We used the MAMBO array at the IRAM 30m telescope to 
study the morphology, masses and densities of the clouds and the dense clumps within them. The latter were followed up by pointed 
ammonia observations with the 100m Efi"elsberg telescope, to determine rotational temperatures and kinematic distances. 
Results. Extinction maps of the Galactic plane trace large scale structures such as the spiral arms. The extinction method probes lower 
column densities, Nu^ ~ 4 x 10"" cm"-, than the 1.2 mm continuum, which reaches up to A'h, ~ 3 x 10'' cm"- but is less sensitive 
to large scale structures. The 1.2 mm emission maps reveal that the high extinction clouds contain extended cold dust emission, from 
filamentary structures to still diffuse clouds. Most of the clouds are dark in 24 /jm, but several show already signs of star formation 
via maser emission or bright infrared sources, suggesting that the high extinction clouds contain a variety of evolutionary stages. The 
observations suggest an evolutionary scheme from dark, cold and diffuse clouds, to clouds with a stronger 1.2 mm peak and to finally 
clouds with many strong 1.2 mm peaks, which are also warmer, more turbulent and already have some star formation signposts. 
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(N 1. Introduction l2006l:lPiUai et alj2006tlSimon et al.l2006bh . From IRDC clumps 

, r , to the next stage, the high-mass protoste llar objects (HMPOs, 

T-H Massive stars play a fundamental role m the evolution of galax- [fieuther et aUM LSridhar an et al.ll2002h . the temperatures in- 

in les through their strong UV radiation, stellar winds and super- ^^^^^^ K < T < 6 K), the hne widths increase, densities and 

^ novae explosions, which contribute to the chemical enrichment ^^^^^^ ^-^^ ( M^Fan[200l . HMPOs are usually found prior 

O of the interstellar medium. Massive stai-s are rare hence usu- formation of ultra compact Hii (UCHii) regions, before the 

O ally found at large distances. They form very rapidly while still ^^^^ ^^^^^^ ^^^^ ^ -^^ -^^-^^ -^^ surrounding medium. 

T— I , deeply embedded in their natal molecular clouds. These char- 

L; ■ acteristics impose several observational obstacles, hke the ne- |Motte et alj (|2007|) demonsti-ated the difficulty of finding 

. ^ ' cessity of high resolution and sensitivity in an un-absorbed fre- massive objects in an early evolutionary phase: in their survey of 

quency range, to study their formation. 3 deg^ in Cygnus X, they found litde evidence for dense clumps 

J_j Currently, the earliest stage of massive star formation is without any trace of star formation, however dense clumps 
, thought to take place in the very dense clumps found in Infrai-ed with already ongoing star formation were found to be abun- 
Dark Cloud s (IRD Cs). The properties of IRDCs are shown by dantly present. Based on these results, the statistical life time 
ICarev etal.1 (if^) to be dense (n^^, > lO'^cm-^) and cool of the high-mass protostars and prestellai- cores was estimated 
(T < 20 K) aggregations of gas and dust in the G alaxy. They con- < 3 x lO"* yr (Motte et al. 2007), which is much shorter than what 
tain clumps with typical masses of > 100 Mq (iRathborne et all is found iri nearby low-mass star-forming regions: ~3 x 10^^ yr 
dKirk et al.l 12005) . It is these very early stages, which provide 



* Tables [T] and [5] are only available in electronic form important clues to construct a theoretical model of massive star 

at http://www.aanda.org Table [T] and fits images associated formation, since the initial fragmentation of the gas and dust in 

with the extinction maps of Fig. □ and the 1.2 mm contin- a clump will be differ ent in the case of monoUthical collapse 

uum maps appearing in Fig. \5\ can be queried from the CDS (Krumholz et al. 2009) co mpared to the com petitive accretion 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via model (B onnell et al.. 200TI:IClark et aljl205l) . 

|http ://cdswed.u-strasbg.fr/cgi-bin/gcat?J/A+A/ 

** Member of the International Max Planck Research School (IMPRS) Massive Stars generally form in clusters (|Lada & Lada| 

for Astronomy and Astrophysics at the Universities of Bonn and 2003), of which the precursors are massive clumps or the so- 

Cologne called precluster forming clumps, hereafter just clumps, of a 
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~ 1 pc size. For many massive star-forming regions, we do not yet 
have the capacity to resolve the clumps into prestellar cores and 
study the fragmentation dBe uther et al. 2007; Rathborne et al] 
120071 l2008t IZhangetal.ll2009t ISwifd 120091) . In this paper, we 
report on the physical parameters of the clumps, such as their 
morphology, density and temperature. Based on this, we hypoth- 
esize on an evolutionary sequence of cluster formation. 

Our understanding of the clumps increased considerably 
with the discovery of IRDCs. IRDCs are detected by a local 
absence of infrared (IR) emission against the diffuse mid IR 
emission of the Galactic plane (Peraultet al. 1996; Ega n et al] 
1 1998) a nd are obser ved numerously throughout the Milky Way 
(Simon et al.ll2006bl) . At the typical low temperatures of IRDCs 
(~ 20 K), the dust emission peaks in the far-infrared and is opti- 
cally thin at mm/submm wavelengths. For a majority of IRDCs 
the mm dust emission coincides wit h the morphology of the 
IR absorption dRathborne etal] 120061; iPillai etal.l l2006V Many 
clumps in IRDCs show signs of star formation via infrared emis- 
sion at 24 um, or S iO emission from shocks dr i ven by outflows 
(Motte e t al J 120071; iBeuflier & Sr idharan' ' 20071; [Chambers et al.1 
[2009). Observations tell that clumps in IRDCs span a very wide 
range of mass es, indicating that not all will form cl usters with 
massive stars dRathborne et al.ll2006tlPiUai et al.ll2006.) . 

The detection method of IRDCs is very sensitive to the lo- 
cal properties of the background emission. Also, not all massive 
dust condensations will be infrared dark if there is enough fore- 
ground emission. Hence, to find the high mass end of molec- 
ular clouds in an unbiased fashion, new, complementary ap- 
proaches are needed. We have developed such a new method, 
well known from studies of low mass star-formi ng regions, 
to targ et more efliciently the most massive clouds: iLada et aTl 
d 19941) pioneered the method of measuring high amounts of ex- 
tinction through stellar color excess in the infrared. Applied to 
the 2/im data of the 2MASS survey, they covered the range 
up to 40 magnitudes in visual extinction. Ay- However, this is 
not sufficient to probe the dense buthplaces of massive stars. 
Here the results of the Spi tzer Space Telescope GLIMPSE sur- 
vey dBeni amin et al."2003') came to help: by applying the ex- 
tinction curve of Indebetouw et al. (2005) we have extended the 
color excess method to reach up to peaks in Ay of ~100 mag- 
nitudes (or column densities of 9 x 10^^ cm"^), thus en- 
tering the realm where massive star formation becomes possi- 
ble. The extinction method, however, is limited by the number 
of available background stars, and will therefore detect mainly 
nearby clouds (discussed in Sect. 12.11 ). In the meanwhile, com- 
plementary, unbiased dust continuum surveys were carried out: 
the ATLASGAL s urvey o f the c omplete inner Galactic plane at 
810 fim by Schuller et al. (2009) and the 1.1mm BOLOCAM 
survey (R osolowsky e t al. 2009) of the Galactic plane accessible 
from the northern hemisphere. 

We selected the more compact and high extinction (mean 
Ay > 20 mag or A^H2 > 2 x 10^^ cm"^) sources from large scale 
extinction maps of the inner Galactic plane (-60° < / < 60°, 
0.9° < b < -0.9°). These high extinction clouds (HECs) were 
studied in the millimeter dust continuum and the rotational tran- 
sitions of ammonia (NH3). NH3 has proven to be a reliable tracer 
of dense gas in dark clouds: not only does the NH3 e missio n 
match the submillimeter dust emission peaks dPiUai et alJl2006l) . 
but also do observations show that, unlike other molecules, NH3 
d oes not deplete fro m the gas phase for typical IRDC densities 
- iPagani et al.l d2005l) observed that NH3 depletes at densities of 
10^ cm "^ in agreement with the prediction of Bergin & Langer 
dl997l) . Moreover, throughout the evolutionary stages of mas- 
sive star formation, NH3 shows an increasing trend in averaged 



line widths and te mperatures from less to more evolved sources 
dPillai et alJ2006h . which indicates that the NH3 molecule is also 
a tracer of evolutionary phase. 

This study presents an overview from the high extinction 
complexes on Galactic size-scales, covering several tens of par- 
sees, to the clumps found in the 1 .2 mm continuum of 0. 1 -0.7 pc 
in size. The connection between the largest and smallest scale is 
important for a comprehensive view of cluster formation in giant 
molecular clouds. In section 2 , we present the method of extinc- 
tion mapping. Observations and data reduction are described in 
section 3, and the results follow in section 4. The analysis and 
discussion of the physical parameters are given in section 5, and 
are compared to previous studies and theoretical predictions in 
section 6. 

2. High extinction clouds 

2.1. Method: extinction mapping 

The overall distribution of dust in a cloud can be traced by the ex- 
tinction of background starlight at visu al and near-infrare d wave- 
lengths as it passes through a cloud dLada et al.lll994l) . Since 
extinction decreases with wavelength, observations at longer 
wavelengths probe deeper into the cloud and trace denser re- 
gions. Additionally the number of detectable background stars 
increases at these wavelengths. With the advance of infrared 
cameras it became possible to detect several hundreds of back- 
ground stars through a cloud, allowing to convert the infrared 
images covering them in extinction maps of useful resolution. 

The GLIMPSE survey employed the Infrared Array Cameras 
dFaz io et al. 2004) onboard the Spitzer Space telescope, operat- 
ing at 3.6, 4.5, 5.8 and 8.0 fim. We used the data provided by the 
GLIMPSE I survey (release April 2005), which covered longi- 
tudes of / = 10 - 65°, 295 - 350° widi b = + 1°. The calibration 
of the data is described in iReach et alJ (l2005h . The inner 20° of 
the Galactic Plane, except for the innermost ±1°, was taken from 
the GLIMPSE II survey (2007). 

To construct the extinction maps we used the averaged 
(3.6;um-4.5;um) color excess, because the extinction law de- 
termination for these wavelengths is the most accurate of all 
the Spitzer bands. The averaged color excess, < £(3.6 ;um - 
4.5 yum) >, was calculated from the color excess in a large scale 
field (a box of size 108" x 108"): 

< £(3.6 fim - 4.5 fim) >- 

< (3.6 i-im - 4.5 fim) - (3.6 //m - 4.5 //m)o > , (1) 

where the background stars are taken to be common-type K gi- 
ants. Measurements of (3.6 fim - 4.5 fim)o in such control fields 
showed that K giants have an average color of ~Omag with 
a dispersion of 0.2 mag. Star ting fr om the reddening law for 

< E(H - K) > in Lada et alj (1 1 994 . one can extend it foflow- 
ing llndebetouw et al.. (.2005) . and get the relation between the 
averaged color excess to the averaged visual extinction, < Ay >: 

< Ay >= 81.8 < £(3.6 /Km -4.5 yum) > [mag]. (2) 

The color excess map can be contaminated by embedded 
stars in the cloud itself or by foreground stars. The latter will 
increase in number as the cloud is located at a farther distance. 
Since the foreground stars will not be reddened, they decrease 
the average color excess of the field. For example, if the num- 
ber of foreground stars equals the number of background stars 
the color excess will be halved. It also means t hat for far away 
clouds the color excess will be underestimated. IRathborne et al.l 
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(120061) and lChambers et all (l2009l) find signs of active massive 
star formation in one of our clouds. Such young red objects will 
contribute to the measured color excess, which will lead to an 
overestimation of the derived extinction. However, the selection 
of clouds associated with very early phases of star formation will 
not be affected. 

In general, the limits of extinction mapping are set by the 
number of available background stars; their number has to be 
sufficient for a statistically meaningful color excess determina- 
tion. Thus, the reach of the extinction method will change with 
Galactic latitude, because at higher latitudes the number of stars 
decreases. In the Galactic plane, there will be "horizon" to which 
one can measure a sufficient color excess, however this hori- 
zon will be far from uniform; it depends for every direction on 
the number of K giants in front and behind the clouds, which 
will differ when crossing a spiral arm or moving in toward the 
Galactic center 

Recently, IChapman et al.l (l2009l) studied the changes in the 
mid-infrared extinction law within a large region with high res- 
olution. They find that while in regions with a K-band extinc- 
tion of Ak < 0.5 mag the extin ction law is wel l fitted by an 
extinction factor of Ry - 3.1 (IWeingartner & Draind 1200 ih . 
the regions with Ak > 1.0 m ag are more consistent with the 
IWeingartner & Draind (1200 Ih model of - 5.5, which uses 
larger maximum dust grain sizes. The high extinction clouds are 
by definition very dense regions for which Ak > 1.0 mag. This 
means that the visual extinctions and column densities reached 
by the extinction mapping are a factor ~1.8 higher than esti- 
mated from Eq. [2] where we used Ry = 3.1. The near and 
mid-infrared extinction law is of great interest, and recently 
many publica tions appeared on the mid-infrared (i Flahertv et al.l 
120071; iNishiv ama et al. 2009; Zasowski et al. 20 0^ and flie 
near-infra red (Moore e t al.i i2005: .Froebrich & del Burgoll200"6t 
IStead & Hoare 2009) extinction law. Within the range of 3.6 and 
4.5 fim, the range we used for our extinction maps, th ere is a 
reasonable agreement between the results of ilndebetouw et al.l 
(I2005h and most recent studies. 

2.2. Catalog of high extinction clouds 

We made extinction maps for the complete inner part of the 
Galaxy with a resolution of 108", plotted in Fig. [1] The mid- 
infrared extinction is changing with longitude and latitude, be- 
cause it is sensitive to large-scale structures. Fig. |2] shows his- 
tograms of the average color excess in longitude and latitude. 
Most of the peaks in the color excess can be associated to 
Galactic spiral arms. The large-scale structure of the extinction 
was excluded by selecting only compact high extinction regions 
with an color excess above 0.25 mag (equivalent to a hydrogen 
column density of 2 x 10^^ cm"^). In a second step, smaller ex- 
tinction maps with a higher resolution (54") were computed to 
obtain more accurate positions of the highest extinction peaks. 
These peaks were selected by eye. Regions with known sources, 
such as Hi: regions and HMPOs, were discarded, leaving a sam- 
ple of unknown and possibly cold and massive clouds. In this 
paper we studied 25 high extinction clouds (HECs) in the first 
Galactic quadrant. These were the clouds, visible from the north- 
ern hemisphere, with the highest extinctio n peaks which were 
not associated with Hii regions from the iBecker et al.l d 19941) 
survey and which had no available mm/submm maps from the 
literature. A complete catalog of all the high extinction clouds 
is given in Table [T] where for each cloud the center position in 
J2000 coordinates is listed with its corresponding peak color ex- 
cess. 



3. Observations and data reduction 

3. 1 . Millimeter bolometer observations and calibration 

The 1.2 mm dust continuum of the selected extinction peaks 
was imaged with the 117-element Max-Planck Millimeter 
Bolometer array (MAMBO-2) installed at the IRAM 30m tele- 
scope. The observations were performed in four sessions be- 
tween 2006 October and 2007 March. The MAMBO passband 
has an equivalent width of approximately 80 GHz centered on 
an effective frequency of 240 GHz {A ~ 1.25 mm). The full 
width at half maximum (FWHM) beam size at this frequency is 
10."5. The maps were taken in the dual-beam on-the-fly map- 
ping mode, where the telescope was scanning row by row in az- 
imuth at a constant speed, while the secondary beam was wob- 
bling in azimuth with a throw of 92". The map sizes were 6'x6', 
resulting in ~ lO'xlO' images with a reduced sensitivity at the 
edges. Each scan was separated by 8" in elevation. We used a 
scanning velocity of 8" s ' with a wobbler rate of 2 Hz. The 
observation time per map was ~ 27 minutes. In total 25 high 
extinction clouds listed in Table |2] were observed, the respective 
center positions used for the observations are given in Table [1] 

The project was observed in the bolometer pool as a backup 
project resulting in strongly varying weather conditions from 
one session to another. The zenith atmospheric opacity, r, was 
monitored by sky dips, performed every ~1.5 hour, and varied 
between 0.1 and 0.5. The calibration was performed every ~3 
hours mainly on H ii regions G34.3H-0.2 and GIO. 2-0.4, and 
was found to be accurate within ~10%. Pointing checks were 
made roughly every half hour to hour, before and usually also af- 
ter each map. Bolometer observations are dominated by the sky 
noise, the variation of the brightness of the sky, which usually 
exceeds the intensity of astronomical sources. The sky noise was 
reduced by subtracting the correlated noise between the bolome- 
ters in the aiTay. The average rm.s. noise signal in the individual 
maps was better than ~15 mJy beam"' after reduction of the sky 
noise. The data were reduced using the MOPSI software pack- 
age developed by R. Zylka. 

3.2. Ammonia observations 

For each cloud, we selected by eye the mm emission peaks for 
follow-up with pointed ammonia observations. All mm emis- 
sion peaks above two times the mean emission of the cloud, 
that is the r.m.s. emission in the cloud where we omit the bright 
mm peaks, were observed. When no clear emission peak was 
present, such as in a very diffuse cloud, the center of the dif- 
fuse emission was targeted. Even in cases where the mm emis- 
sion was below 3cr, where cr is the noise in the bolometer map 
(cr - 0.015 Jy beam"'), we choose a few positions to observe 
the ammonia lines to search for a cloud so cold and diffuse that 
is was missed by the bolometer For clouds that had several con- 
tinuum emission peaks separated by more than 20" from each 
other, more than one position was observed, denoted by MMl, 
MM2 etc. Several high exinction clouds had very weak or no 
mm emission - in this case the name extension, which was MMl, 
MM2 for the mm sources, was changed to 1 or 2, e.g., G034.34- 
00.90 1. 

The ammonia observations were performed with the MPIfR 
100m Effelsberg Telescope. We observed 54 positions, given in 
Table [3] and g] between 2007 April and 2008 December. The 
Effelsberg beam is 40" FWHM at the ammonia inversion line 
frequencies of ~23.7 GHz. The 2007 observations used the 8192 
channel AK90 auto correlator backend. The correlator was con- 
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figured into eight spectral windows with 20 MHz bandwidth 
and 1024 channels each, where every window could be set to 
a different frequency. This provided the opportunity to simul- 
taneously observe the {J,K) = (1,1), (2,2) and (3,3) inversion 
lines of ammonia in both polarizations. The spectral resolution 
with this setup was 0.25 km s"' . The 2008 observations used the 
Fast Fourier Transform Spectrometer (FFTS) with a 500 MHz 
bandwidth and 8192 channels. This bandwidth is sufficient to si- 
multaneously observe the three ammonia inversion lines with a 
spectral resolution of 0.7 km s"'. All the observations were per- 
formed in frequency switching mode with a throw of 7.5 MHz. 



3.2.1. Data calibration 

During the observations, pointings on a nearby compact contin- 
uum source were performed every hour for determining point- 
ing corrections. For the flux calibration, we observed a well- 
known flux calibrator, NGC 7027 or 3C 286, in every run. The 
post-observational calibration to obtain the main beam bright- 



ness temperature. Tub, consisted of the opacity correction, ele- 
vation correction, and flux calibration: 



' MB 



Giff) 



(3) 



where A is a scaling factor, is the antenna temperature, r the 
zenith opacity, and G the function of the gain with elevation 0. 

The opacity, t, was calculated by fitting a linear function to 
the system temperature against airmass (sin"' 6) and taking the 
slope of the fit. We found a r of 0.03 1 for good weather When 
it was not possible to retrieve the r by this method the average 
T of 0.054 at 23 GHz was assumed (the averaged water vapor 
radiometer value for 2007). 

For all parabolic dish telescopes the gain decreases at very 
low and very high elevations. We corrected for this by dividing 
by G(ff) taken from the Effelsberg websit^ (see Eq.[3]), which is 
given by 



G(0) - a() + aiO + a2(r 



(4) 



http://www.mpifr-bonn.mpg.de/div/effelsberg/calibration/l. 3cmsf.html 
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HEC name 




Radius " 


< Av >cxt 


< A'h, >exl 


F' 




Jyl 1 T nun 

1 .z i[]in 


Clump/Cloud^ 


Class 






("/pc)~ 


(mag) 


(10^2 cm-2) 


(Jy) 


(Mo) 


(Mo) 


(mjy beam^'/rnJy) 




G012.73-00.58... 




109/0.26 


32 


3.0 


1.89 


119 


109 


63/52 


diffuse 


G013.28-00.34... 




59/1.15 


42 


3.9 


3.85 


3,039 


1,693 


129/66 


diffuse 


G013.91-00.51... 




49/0.65 


44 


4.2 


2.84 


1,015 


602 


203/66 


peaked 


G013.97-00.45... 




52/0.6 


22 


2.1 


2.29 


457 


308 


99/56 


diffuse 


G014.39~00.75... 


A 


62/0.66 


34 


3.2 


4.88 


507 


796 


64/195 


peaked 




B 


25/0.31 


16 


1.5 


0.69 


86 


165 


103/66 


diffuse 


G014.63-00.57... 




98/1.04 


35 


3.3 


17.47 


2,084 


1,822 


907/96 


multiply p. 


G016.93+00.24... 




52/0.61 


32 


3.0 


2.20 


632 


294 


118/56 


peaked 


G017. 19+00.81... 




71/0.79 


34 


3.2 


6.78 


1,258 


785 


547/63 


multiply p. 


G018.26-00.24... 




85/1.93 


32 


3.0 


11. n 


6,540 


5.896 


317/92 


multiply p. 


G022.06+00.21... 




45/0.81 


28 


2.6 


4.05 


1,000 


1,018 


1151/72 


multiply p. 


G023. 38-00. 12... 




68/1.84 


29 


2.8 


5.16 


5,498 


3,470 


211/67 


peaked 


G024.37-00.15... 




43/0.81 


26 


2.5 


0.56 


921 


561 


168/52 


multiply p. 


G024.61-00.33... 




40/0.61 


25 


2.4 


1.78 


496 


424 


214/69 


multiply p. 


G024.94-00.15... 




49/0.79 


32 


3.0 


2.17 


1,095 


650 


160/62 


peaked 


G025.79+00.81... 




38/0.63 


24 


2.3 


1.05 


494 


351 


-/46 


diffuse 


G030.90+00.00 ... 


A 


36/0.80 


24 


2.2 


1.45 


826 


619 


188/61 


peaked 




B 


39/1.35 


20 


1.8 


1.57 


1,960 


2,232 


157/57 


peaked 




C 


30/0.83 


21 


2.0 


0.83 


799 


636 


129/71 


diffuse 




D 


14/0.17 


22 


2.1 


0.17 


37 


31 


115/68 


diffuse 


G034.03~00.33... 




6/- 


21 


1.9 


0.02 


- 


- 


-/60 


diffuse 


G034.34~00.90... 




-/- 


- 


- 


- 


- 


- 


-/- 


diffuse 


G034.71~00.63... 




51/0.74 


31 


3.0 


3.05 


968 


642 


206/54 


multiply p. 


OL/jH". / /— UU.Oi... 




Q/n 1 




^ n 

J.U 


U.UJ 


j\j 


1 9 


_j_ 


diffuse 


G034.85+00.43... 




16/0.29 


37 


3.5 


0.17 


171 


73 


-/43 


diffuse 


G035.49~00.30... 


A 


19/0.34 


11 


1.1 


0.39 


73 


103 


198/58 


peaked 




B 


63/0.91 


33 


3.1 


4.1 


1,495 


1,300 


118/56 


peaked 


G037.44+00.14... 


A 


26/0.16 


33 


3.1 


0.58 


47 


37 


-/67 


diffuse 




B 


9/0.11 


25 


2.4 


0.05 


17 


9 


-/54 


diffuse 


G050.06+00.06... 




51/1.18 


29 


2.8 


2.92 


2,257 


2,227 


140/54 


peaked 


G053.81~00.00... 




40/0.37 


31 


2.9 


1.54 


236 


192 


140/56 


peaked 



Notes. The first two columns give the HEC name and the subgrouping of clouds based on the kinematic distance, the following columns represent 
(in order of appearance): radius of the cloud, average visual extinction derived from the extinction maps, average hydrogen column density derived 
from the extinction maps, integrated 1.2 mm flux, cloud mass derived from the extinction maps, cloud mass derived from the 1.2 mm, clump peak 
flux and median cloud flux, class of HEC. 

" Radius is the square root of the area divided by n. For the conversion to parsecs we use the kinematic distance based on the NH3(1,1) line. 
* Peak of the 1.2 mm emission in the clump and the mean 1.2 mm emission of the cloud without the clumps. 



where oq = 0.88196, oi = 6.6278 x 10"^ 02 = -9.2334 x 
10"^, and 9 the elevation. After the opacity and gain correc- 
tions, the scaling factor A was found by comparing the mea- 
sured uncalibrated intensities of absolute flux density calibrator 
sources with the lite r ature values calculate d from the form ulae 
given bv lBaars et all (Il977h for 3C 286, and lOtt et alJ (11994 for 
NGC 7027. 



3.3. Water maser observations 



The calibrated spectra were baseline subtracted, and the am- 
monia lines were fitted by a Gaussian. Only the NH3(1,1) line, 
for which the hyperfine structure was clearly detectable given 
the signal to noise ratios, was fitted by special routine 'method 
nh3(l,l)' of the GILDAS/CLASS software. This method calcu- 
lates the optical depth from the hyperfine structure and returns 
optical depth coiTected line widths. The observed ammonia pa- 
rameters, such as the velocity in the local standard of rest (LSR), 
Vlsr^ main beam temperatures, Tmb, line widths, Av, and the 
main group optical depth, T^ain are listed in Table|5j 



Several positions with a peak in the 1.2 mm emission above 
twice the mean cloud emission, indicating that they are possi- 
bly harboring evolved clumps, were searched for water maser 
emission using the 100m Effelsberg telescope on 5 and 24 of 
February 2008. In total 24 positions in common with NH3 were 
observed (marked in Table[3]). We performed on-off observations 
using the FFTS backend with a bandwidth of 20 MHz centered 
on 22.235 GHz. This setup afforded a high spectral resolution of 
0.04km s"', while allowing the water maser to have line widths 
of up to 100km s"'. The bandwidth was 270km s"' and was 
centered on the NH3 Vlsr of the observed clump. We applied 
the same data reduction as described above in Sect. 3.2.1. for 
the ammonia observations. Table |6] lists the Vlsr and peak in- 
tensity of the water maser detections. In case of multiple maser 
components we give the Vlsr range. 
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Table 3. Positions (J2000) and fluxes of the clumps from the 1.2 mm continuum data 



HEC name 




R.A. 


Declination 


pr 




(mJy) 


Major Axis 


Minor Axis 


P.A. 


Hzo;;, 






(h:m;s) 


(°:':") 


(mJy beam" ^ ) 


(mJy) 


(arcsec) 


(arcsec) 


(deg) 




G012.73-00.58... 


MMl 


18:15:41.3 


-18:12:44 


42(8) 


165 


316 


28.4 


15.2 


26.9 




G013. 28-00.34... 


MMl 


18:15:39.9 


-17:34:37 


59(12) 


212 


141 


22.7 


17.5 


-49.6 


- 


G013.91-00.51... 


MMl 


18:17:34.8 


-17:06:52 


142(15) 


685 


535 


27.3 


17.1 


-22.2 


- 


G014.39-00.75A.. 


MMl 


18:19:19.0 


-16:43:49 


142(17) 


503 


456 


27.1 


14.4 


57.1 


- 


G014.39-00.75B.. 


MM3 


18:19:33.3 


-16:45:01 


81(9) 


259 


2,11 


21.5 


16.4 


6.7 




G014.63-00.57... 


MMl 


18:19:15.2 


-16:29:59 


689(8) 


3,111 


2.780 


28.5 


17.5 


10.3 


+ 




MM2 


18:19:14.3 


-16:30:41 


332(9) 


584 


650 


16.2 


12.0 


-57.3 






MM3 


18:19:02.9 


-16:30:29 


112(11) 


269 


276 


21.5 


12.3 


-42.0 






MM4 


18:19:20.5 


-16:31:42 


86(9) 


215 


265 


22.1 


12.5 


-49.3 




G016.93+00.24... 


MMl 


18:20:50.8 


-14:06:01 


76(11) 


357 


296 


27.2 


19.1 


-80.4 




G017. 19+00.81... 


MMl 


18:19:08.9 


-13:36:29 


120(10) 


386 


409 


21.6 


16.4 


65.3 






MM2 


18:19:12.9 


-13:33:46 


475(11) 


1,304 


1,430 


17.8 


16.9 


87.8 


+ 




MM3 


18:19:12.1 


-13:33:32 


157(11) 


725 


1,080 


31.3 


16.2 


-72.9 






MM4 


18:19:15.2 


-13:39:29 


87(10) 


426 


407 


28.2 


19.1 


-42.0 




G018.26-00.24... 


MMl 


18:25:11.8 


-13:08:04 


257(16) 


666 


438 


19.4 


14.7 


26.5 


- 




MM2 


18:25:06.4 


-13:08:51 


221(17) 


887 


429 


31.2 


14.2 


30.2 


- 




MM3 


18:25:05.6 


-13:08:20 


152(16) 


231 


201 


14.8 


11.3 


-42.8 






MM4 


18:25:04.5 


-13:08:27 


133(14) 


211 


168 


16.7 


10.5 


47.7 






MM5 


18:25:01.8 


-13:09:06 


141(19) 


559 


291 


24.1 


18.1 


64.6 




G022.06+00.21... 


MMl 


18:30:34.7 


-9:34:46 


995(14) 


2,019 


1,820 


15.6 


14.3 


70.6 


+ 




MM2 


18:30:38.5 


-9:34:29 


166(15) 


259 


208 


14.6 


11.8 


64.0 




G023.38-00.12... 


MMl 


18:34:23.5 


-8:32:20 


171(15) 


521 


246 


22.3 


15.1 


29.1 




G024.37-00.15... 


MMl 


18:36:27.8 


-7:40:24 


133(12) 


529 


350 


26.8 


14.4 


33.0 


- 




MM2 


18:36:18.3 


-7:41:00 


94(11) 


254 


202 


17.3 


15.1 


42.7 


+ 


G024.61-00.33... 


MMl 


18:37:23.1 


-7:31:39 


147(8) 


655 


496 


25.3 


19.4 


84.2 


- 




MM2 


18:37:21.3 


-7:33:07 


75(8) 


276 


187 


26.1 


15.6 


20.6 


- 


G024.94-00.15... 


MMl 


18:37:19.7 


-7:11:41 


133(15) 


350 


318 


19.4 


13.2 


-32.8 


-1- 




MM2 


18:37:12.2 


-7:11:23 


114(15) 


278 


218 


19.2 


12.4 


23.2 


- 


G030.90+00.00A.. 


MMl 


18:47:28.9 


-1:48:07 


141(10) 


728 


307 


32.1 


17.7 


-2.5 




G030.90+00.00B.. 


MM2 


18:47:41.9 


-1:52:13 


116(19) 


494 


130 


26.9 


17.4 


43.7 


- 


G030.90+00.00C.. 


MM3 


18:47:48.2 


-1:51:30 


98(20) 


369 


164 


24.3 


17.1 


-10.5 


- 


G030.90+00.00D.. 


MM4 


18:47:51.5 


-1:49:24 


106(10) 


260 


279 


20.9 


13.0 


32.4 




G034.71-00.63... 


MMl 


18:56:48.3 


1:18:49 


148(9) 


1,174 


678 


36.0 


24.3 


2.2 


- 




MM2 


18:56:58.2 


1:18:44 


79(11) 


368 


250 


30.2 


17.1 


-81.7 






MM3 


18:57:06.5 


1:16:52 


60(11) 


203 


179 


23.9 


15.6 


78.9 


- 


G035.49-00.30A.. 


MMl 


18:57:05.2 


2:06:29 


180(11) 


374 


339 


16.8 


13.6 


-38.8 


+ 




iVliVlZ 


lO.D 1 .Uo.H- 








1 


1 Q 1 




Zo.U 






MM3 


18:57:08.1 


2:10:47 


73(14) 


437 


279 


33.3 


19.9 


-1.6 


-1- 




MM4 


18:57:09.0 


2:08:23 


60(10) 


178 


118 


24.1 


13.5 


-1.4 






MM5 


18:57:06.7 


2:08:27 


62(10) 


211 


154 


26.4 


14.3 


-28.9 






MM6 


18:57:11.5 


2:07:27 


55(11) 


274 


184 


29.0 


19.1 


-6.1 




G050.06+00.06... 


MMl 


19:23:12.4 


15:13:35 


89(8) 


332 


167 


26.9 


15.4 


43.1 






MM2 


19:23:09.2 


15:12:42 


84(9) 


246 


157 


30.5 


14.8 


83.0 




G053.81-00.00... 


MMl 


19:30:55.7 


18:29:55 


106(11) 


226 


291 


16.5 


14.3 


-50.8 





Notes. The first two columns give the HEC name and the millimeter clump number, the following columns represent (in order of appearance): 
right ascension, declination, 1.2 mm peak flux, integrated 1.2 mm flux, integrated 1.2 mm flux within 0.25 pc diameter, clump major axis, clump 
minor axis, position angle and water maser detection. 



" Water maser; '+'means a detection, '-' a non-detection and not observed. 



4. Results 

4. 1 . Kinematic distances 

Accurate distance determination within the Galaxy is generally 
difficult. Nevertheless, the kinematic distance is commonly used 
as a distance measure. It is based on a model of Galactic rota- 
tion, the "rotation curve" characterized by the distance between 
the Sun and the Galactic Center, Rq, and the rotation velocity at 
the Suns orbit, ©q. Toward the inner part of the Galactic plane, 
kinematic distances are ambiguous: for a given Galactic longi- 
tude and LSR velocity, it cannot a priori be determined if the 
object is at the "near" or "far" kinematic distance. All extinc- 



tion clouds should be at the near kinematic distance since at 
large distances the number of background stars decreases and 
the percentage of foreground stars increases making it difficult to 
measure any color excess. We calculated the kinematic distances 
for all clumps with NH3 detections, using a progra m of Todd 
Hunter, which applies the Galactic rotation model of iFich et al.l 
(1989) assuming a flat rotation curve, ©o = 220 kms"', and 
Rq = 8.5 kpc. The resulting kinematic distances are given in 
Table|4]and|2l 
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Table 4. Positions (J2000) and determined properties toward the positions without clumps 



7 



HEC name 






R.A. 


Declination 




-^^0 25 


^0.25 pc 








(h:m:s) 


(":':") 


(kpc) 


(mJy) 


(Mo) 


(K) (10'^ cm-^) 


G012.73-00.58... 


MM2 


18 


15:32.7 


-18:10:15 


1.1 


178 


8 


11.4(0.7) 


G013.97~00.45... 


MMl 


18 


17:16.5 


-17:01:16 


2.4 


189 


25 


16.6(0.9) 


G014.39-00.75A.. 


MM2 


18 


19:17.4 


-16:44:04 


2.1 


382 


32 


20.0(2.3) 


G023.38-00.12... 


MM2 


18 


34:20.4 


-8:33:16 


5.6 




44 


17.8(2.0). 


G025.79+00.81... 


MMl 


18 


35:20.5 


-5:56:36 


3.4 


196 


77 


13.1(0.7) 




iVliVlZ 


18 


35:26.3 






Oj 


^1 

J i 




G034.03-00.33... 


MMl 


18 


54:25.1 


0:49:56 












2 


18 


54:39.2 


0:51:37 










G034.34-00.90... 


1 


18 


57:16.5 


0:50:48 










G034.77~00.81... 


MMl 


18 


57:40.7 


1:16:09 


2.9 


49 


11 




G034.85+00.43... 


MMl 


18 


53:23.2 


1:53:16 


3.6 


73 


31 


13.2(1.5) 


G037.44+00.14A.. 


MMl 


18 


59:14.0 


4:07:37 


1.3 


163 


8 




G037.44+00.14B.. 


MM2 


18 


59:10.2 


4:04:32 


1.7 


130 







Notes. The first two columns give the HEC name and the millimeter clump number, the following columns represent (in order of appearance): 
right ascension, declination, kinematic distance, integrated 1.2 mm flux within 0.25 pc diameter, mass within 0.25 pc diameter, NH3 rotational 
temperature, and NH3 column density. 



Table 6. Detected water masers 



HEC name 


Peak Intensity 


Vlsr 




(Jy beam"') 


(km s-') 



G014.63-00.57 


MMl 


9.3 


22 


G017.19+00.81 


MM2 


22.3 


-4-37 


G022.06+00.21 


MMl 


11.7 


42-51 


G024.37-00.15 


MM2 


6.0 


65 


G024.94-00.15 


MMl 


6.6 


52-72 


G035.49-00.30A 


MMl 


5.6 


40-77 


G035.49-00.30B 


MM2 


2.6 


40 



4.2. Galactic scale extinction 

The low resolution extinction maps show large extinction struc- 
tures across the whole inner Galactic disk from 60° to -60° 
longitude (Fig. [TJ. In longitude, the 3° averaged distribution 
(Fig. |2] top panel) shows signs of Galactic structure, meaning 
that peaks in the distribution can be related to known spiral 
arms. The column density appears higher in the fourth quadrant, 
0° > / > -60°, than in the first, 60° > / > 0°. Additionally, there 
seems to be a quasi symmetrical distribution around the Galactic 
Center around longitudes of -15°, -10°and -35°, -30°. In lati- 
tude (Fig. |2] lower panel), we found a peak towards b 0.°1, 

which is also seen for the compact submil limeter sources found 
in the ATLASGAL survey (ISchuUer et al.l20 09). The shift of the 
peak out of the midplane indicates that the Sun is located above 
the midplane. Similar results have been found b y the stu dies us- 
ing young open clusters and OB stars (see e.g. Joshi 2007). The 
FWHM of our distribution is -1°, while in ATLASGAL this is 
more narrow, ~0.°6. If one assumes that the Galactic disk has 
a constant scale height, then objects closer to the Sun should 
have a wider latitude distribution then ones at larger distances. It 
therefore appears that the high extinction clouds are on average 
closer to the Sun than the ATLASGAL sources. 

With the kinematic distances we can place our sample of 
clouds within a face-on view of the Galactic plane (Fig.|3]l. They 
inh abit similar re gions as IRDCs (see the IRDCs distribution 
by iJackson et al.l 12008). The extinction method misses nearby 
(d < 1 kpc) and most of the far away (d > 4 kpc) clouds. The 




Gulactic Latitude [degree] 



Fig. 2. Mean color excess along Galactic longitude (upper dia- 
gram, bins of 3°) and Galactic latitude (lower diagram, bins of 
0.1°). In the longitude histogram, encounters with a spiral arm 
are marked with letters representing N for Norma-Cygnus, Sc 
for Scutum-Crux, Sa for Sagittarius-Carina, and P for Perseus. 



insensitivity to the latter stems from the increasing number of 
foreground stars at larger distances (see Sect. 12. lb . Between dis- 
tances of one and four kilo parsec the high extinction clouds 
agree with the IRDCs regions. 

4.3. Extinction clouds 

Millimeter continuum emission was detected in 88% of the 
bolometer maps. It usually followed well the 54" resolution ex- 
tinction maps (Fig.|4]i. Two examples of weak millimeter emis- 
sion or non-detections are G034.03-00.33 and G034.34-00.90 
(Fig- EJ. Figure |5] shows the 24yum maps from the Multiband 
Imaging Photometer for Spitzer (MIPS) overlaid with the mm 
emission. While the high extinction clouds match well with 
24yum dark features, the 1.2 mm peaks often contain weak com- 
pact lAfim emission. The bolometer maps were sensitive to 
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HEC GO 13. 28-00, 34 



.S-iB°i2'oa" 




R.A. (JHOOO) 



HEC G014-39-00.75 



18 17 20 15 
R.A. (JSOOO) 



-16°44'00" 



-M 



-16°45'00" 



-16°46'00 

16''19"'35^ 




MJy/ sr 
500 



400 



300 



- 200 



100 



30 



25 20 
R.A. (J3000) 



15° 



Fig.5. Spitzer/MIPS 24 //m images overlaid with the 1.2 mm continuum emission contours starting at 2 V2cr (cr - O.OlSJybeam ') 
and increasing by factors of V2. The 1.2 mm clumps are labeled with their designation (as listed in Table [3] Q and |5]). In the bottom 
left corner we indicate the original MAMBO beam at 1.2mm (10."5), or the smoothed beam (20") when a smoothed map is shown. 
Additionally, 6.7 GHz methanol masers (cyan triangles), water masers (blue squares), Hii regions (green pentagons with black 
contours), and IRAS sources are marked in the images. The 1.2 mm continuum emission maps are available in fits format on the 
CDS. 
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HEC G014. 63-00-57 



MJy/sr 
500 



IIRAS 18164-1632 
-\ MMl , 

Fm 14-63-0. 59\r 



IRAS 18164-1631 



— -16°3r00" 



o-16°32'00" 




MM2 



« 



o 



HEC G016, 93+00-24 



10 

R.A. (J2000) 

MJy/ sr 
300 



HEC G017-19 + 00-ai 




-13°10'00' 



R.A. (J2000) 



Fig. 5. -Continued. 



structures from 10."5 (FWHM) to ~ 90"; large-scale structures We find single clouds and clusters of clouds, for example 
are not faithfully represented due to the sky noise subtraction G014. 63-00. 57 in Fig. |5] Based on the kinematic distance the 
and chopping. We smoothed the maps, by convolving with a clouds were identified as belonging together in one large cloud, 
20" Gaussian beam, to study the extended cloud structure and or being physically separated clouds if the difference in distance 
increase the signal-to-noise toward weak and diffuse sources. was larger than 10%. The latter clouds were marked as A, B etc. 

and were treated as distinct clouds throughout the paper. For the 
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Fig. 5. -Continued. 



18 37 ao 15 
R.A. (J200Q) 



total integrated flux of the cloud (Table |2]l, we defined a cloud 
edge at a threshold of 3o- (of the smoothed map) and integrated 
all the emission within this region using the Gildas Software 
package GREG. The corresponding radius of the cloud, VA/tt 
with A the surface of the cloud, is given in arcseconds and par- 
sees (using the kinematic distance) in Table|2] The average cloud 



radius was 0.7 pc. Certain clouds have spherical shapes, for ex- 
ample G013. 28-00. 34, however, most are filamentary, as can be 
seen in G024.94-00.15 or G016.93+00.24. Cloud sizes are very 
diverse depending on the geometry; there are elongated struc- 
tures up to several parsecs in length. The number of clumps 
per cloud varies from five (in G018. 26-00. 24, Fig.|5]l to one or 
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Fig. 5. -Continued. 



none in very diffuse clouds (G034.85+00.43, Fig. |5]l. We used 4.4. Clumps in high extinction clouds 
the SIMBAD database to investigate if our clouds contain signs 

of star formation, such as maser emission, Hn regions, and IRAS ^^^^ ^^^^ unsmoothed maps, as the ex- 

sources 

tended structure suiTounding the clump complicates the defini- 
tion of source edge and confuses Gaussian fitting routines. We 
removed the extended emission from the map by applying a 
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HEC G034.77-00,B1 



EEC G034,B5 + 00.43 




K.A. (JHOOO) 



Fig. 5. -Continued. 



median filter with a box size of 63" (6 x beamsize) using the 
MIRIAD task 'immedian' . After median removal, the source- 
find algorithm of MIRIAD, 'sfind', was run on these images de- 
livering 2D Gaussian fits, the peak flux, and the integrated fluxes 
of the clumps. The median removal should strip the low density 
mass reservoir which surrounds the denser part of the clump. 



First, we selected the peaks in the mm emission by eye, and 
then we took the parameters for the FWHM size, the integrated 
flux and peak flux of the clump from 'sfind'. The observed pa- 
rameters of the clumps are presented in Table [3] The clumps 
have sizes (geometrical mean of the major and minor axis of the 
clump) of 10 - 30" or 0.11 - 0.72 pc, peaking around 0.25 pc 
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Fig. 7. Distance as a function of clump mass (left) and size (right). In the left panel, the clump mass based on Gaussian fitting 
is marked by squares, and the clump mass within 0.25 pc around the center of a clump is indicated by crosses. There is a linear 
correlation with distance for the mass and the size. The most massive clumps are at large distances, whereas almost no low-mass 
clumps are found at distances larger than 4 kpc, which can be expected for a sensitivity limited selection of clouds. 



o 

Q_ 




2 4 

X [kpc] 

Fig. 3. A face-on view of the Galactic plane. The Galactic cen- 
ter is located at (0,0). The Galacitic longitudes of 12° and 53° 
are indicated by solid lines. The locations of the sample of high 
extinction clouds, which were studied in this paper, are marked 
by open circles. 



(Fig. |6]l. Most of the clumps are resolved with the 10."5 beam 
of the 30m telescope. The (physical) clump size follows a linear 
trend with distance (Fig.|7] right panel). At far distances we ob- 
serve only larger, hence brighter clumps, which is expected for 
our sensitivity limited selection. The angular sizes, however, are 
more equally distributed (see Table [3]). For several clouds with 
very weak mm emission the selected position was not the mm 
emission peak, but the emission center. In such clouds 'sfind' did 
not find any clumps; these positions without clumps are listed in 
Table m 



Ammonia emission was detected toward 94% of the ob- 
served positions. Figure [8] shows the spectra of the NH3(1,1), 
NH3(2,2) and NH3(3,3) transitions toward two compact clumps 
representing typical values for an early stage where the clump is 
cold (14K, G016.93-i-00.24 MMl) and a more evolved stage, 
where the temperature has increased (19 K, G017. 19-1-00.81 
MM2). In the early stage the NH3(3,3) line is not detected, while 
for the warmer clump it is clearly present. The results of the 
ammonia observations are summarized in Table |5] The main 
beam brightness temperatures of the NH3(1,1) lines are, after 
baseline removal, between 0.5 and 4.0 K. The baseline r.m.s. 
is ~0.2K, but several spectra are more noisy. The line widths 
from the main component of the NH3 (1,1) range between 0.7 and 
2.8 kms"', yielding an average of 1.4kms"'(Fig.|9l). These line 
widths are far above a thermal linewidth, which would be around 
0.23 km s ' for temperatures of ~18 K. For most of the sources, 
the NH3(1,1) and NH3(2,2) lines are both detected, while the 
NH3(3,3) line is often very weak or not present. The NH3(2,2) 
and NH3(3,3) lines are on average both wider than the main 
NH3(1,1) line. This implies that these lines do not trace exactly 
the same volume of gas, meaning that the beam filling factor is 
not identical. 



5. Analysis 

5.1. Extinction masses 

Since the color excess, or extinction, is a direct measure of 
the amount of column density in a region, we could derive ex- 
tinction masses. One magnitude of Ay is related to a hydro- 
gen column density, A^h, , according to iBohlin et al.i (J 978.) and 
iFrerking et all (1 19821) by 



0.94 X 1021 



[mag]. 



(5) 



Multiplying by the mass of a hydrogen atom, nin, the mean 
molecular weight, fi = 2.33, and the cloud surface. A, one ar- 
rives at the extinction mass: 

Mext = 3.7xl0^^fimu < £(3.6//m-4.5/im) > (^) [Mo].(6) 

^pc"^' 
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Table 7. Properties of the clumps derived from NH3 and the 1 .2 mm continuum. 



HEC name 




,/ 

"kin 


Ay 


r W MM 






-'•^1.2 mm 






T 

rot 








(kpc) 


(mag) 


(pc) 


1 1 n22 2\ 

(lU cm ) 


(iU cm ) 


(Mo) 


(Mo) 


(Mo) 


\ 

(K.) 


(iU cm ) 


G012.73-00.58.. 


MMl 


1.1 


47 


0.11 


4.5 


2.8 


12 


22 


6 


9.3(0.7) 


3.7(0.8) 


G013. 28-00.34.. 


MMl 


4.0 


31 


0.39 


2.9 


0.5 


93 


62 


104 


14.9(2.0) 


3.2(1.1) 


G013.91-00.51.. 


MMl 


2.7 


82 


0.28 


7.7 


2.2 


145 


113 


50 


14.1(1.3) 


2.5(0.5) 


G014.39-00.75A.. 


MMl 


2.2 


61 


0.21 


5.7 


2.0 


52 


47 


31 


17.4(2.0) 


1.1(0.6) 


G014.39-00.75B.. 


MM3 


2.5 


61 


0.23 


5.7 


1.8 


62 


90 


19 


11.8(1.5) 


2.4(0.8) 


G014.63-00.57.. 


MMl 


2.2 


280 


0.24 


26.3 


in 


323 


289 


82 


18.1(1.3) 


4.7(0.5) 




MM2 


2.2 


164 


0.15 


15.4 


7.3 


72 


80 


26 


15.7(1.1) 


2.7(0.4) 




MM3 


2.1 


55 


0.17 


5.2 


2.2 


31 


32 


34 


15.8(1.9) 


3.4(0.8) 




MM4 


2.3 


32 


0.18 


3.0 


1.2 


21 


26 


12 


19.1(5.9) 


2.0(1.1) 


G016.93+00.24.. 


MMl 


2.4 


44 


0.27 


4.2 


1.1 


64 


53 


23 


14.0(1.3) 


1.7(0.3) 


G017. 19+00.81.. 


MMl 


2.5 


52 


0.23 


4.9 


1.5 


54 


58 


35 


17.2(0.9) 


2.1(0.2) 




MM2 


2.3 


186 


0.20 


17.4 


6.3 


142 


156 


35 


18.7(1.1) 


2.1(0.2) 




MM3 


2.3 


56 


0.25 


5.3 


1.5 


72 


108 


45 


20.0(1.5) 


2.6(0.3) 




MM4 


2.3 


31 


0.26 


2.9 


0.8 


40 


39 


211 


20.1(2.1) 


3.2(1.1) 


G018.26-00.24.. 


MMl 


4.7 


104 


0.39 


9.8 


1.8 


306 


201 


91 


18.2(1.1) 


4.3(0.4) 




MM2 


4.7 


95 


0.48 


8.9 


1.3 


431 


208 


201 


17.4(1.2) 


6.5(0.6) 




MM3 


4.7 


76 


0.29 


7.1 


1.7 


130 


113 


150 


15.7(1.5) 


6.4(0.9) 




MM4 


4.7 


61 


0.30 


5.7 


1.3 


110 


88 


140 


16.6(1.1) 


5.3(0.6) 




MM5 


4.6 


63 


0.47 


6.0 


0.9 


278 


145 


197 


16.8(1.7) 


6.0(1.0) 


G022.06+00.21.. 


MMl 


3.7 


270 


0.27 


25.4 


6.8 


384 


346 


81 


24.7(3.3) 


4.9(0.9) 




MM2 


3.7 


84 


0.24 


7.9 


2.4 


92 


74 


36 


15.5(1.6) 


2.6(0.5) 


G023. 38-00. 12.. 


MMl 


5.6 


70 


0.50 


6.5 


0.9 


346 


163 


134 


18.1(1.2) 


4.4(0.5) 


G024.37-00.15.. 


MMl 


3.9 


52 


0.38 


4.9 


1.1 


165 


110 


174 


18.6(1.9) 


7.1(1.2) 




MM2 


3.8 


47 


0.30 


4.4 


1.2 


94 


74 


70 


15.7(1.5) 


3.8(0.7) 


G024.61-00.3.3. 


MMl 


3.1 


63 


0.34 


5.9 


1.2 


140 


106 


69 


17.5(1.9) 


2.0(0.4) 




MM2 


3.2 


38 


0.31 


3.5 


0.8 


72 


49 


21 


15.5(1.2) 


1.5(0.2) 


G024.94-00.15.. 


MMl 


3.3 


69 


0.26 


6.5 


2.0 


104 


94 


61 


15.2(1.0) 


3.8(0.4) 




MM2 


3.4 


59 


0.25 


5.5 


1.8 


85 


67 


52 


15.2(1.1) 


3.5(0.5) 


G030.90+00.00A.. 


MMl 


4.6 


55 


0.53 


5.2 


0.7 


311 


131 


125 


18.6(1.7) 


3.1(0.6) 


G030.90+00.00B.. 


MM2 


7.2 


61 


0.76 


5.8 


0.5 


702 


736 


497 






G030.90+00.00C.. 


MM3 


5.7 


46 


0.56 


4.3 


0.5 


283 


126 


71 


16.4(2.4) 


3.4(0.9) 


G030.90+00.00D.. 


MM4 


2.6 


101 


0.21 


9.5 


3.2 


49 


95 


22 






G034.71-00.63.. 


MMl 


3.0 


61 


0.43 


5.8 


1.0 


224 


129 


130 


17.8(1.1) 


2.2(0.3) 




MM2 


3.0 


55 


0.33 


5.2 


1.1 


122 


83 


59 


12.4(0.9) 


4.4(0.7) 




MM3 


3.1 


26 


0.29 


2.5 


0.6 


43 


38 


59 


17.1(1.9) 


1.0(0.4) 


G035.49-00.30A.. 


MMl 


3.6 


71 


0.26 


6.6 


1.8 


98 


89 


100 


18.6(2.9) 


1.7(1.0) 


G035.49-00.30B.. 


MM2 


3.0 


57 


0.25 


5.3 


1.5 


72 


49 


26 


11.9(0.7) 


3.9(0.5) 




MM3 


3.0 


44 


0.38 


4.2 


0.8 


127 


81 


102 


13.6(0.7) 


4.3(0.4) 




MM4 


3.0 


39 


0.26 


3.7 


1.0 


54 


36 










MM5 


3.0 


40 


0.28 


3.8 


1.0 


64 


46 










MM6 


3.0 


36 


0.34 


3.3 


0.7 


83 


55 








G050.06+00.06.. 


MMl 


4.8 


47 


0.47 


4.4 


0.7 


208 


104 


83 


14.9(2.2) 


2.1(0.6) 




MM2 


4.8 


49 


0.49 


4.6 


0.7 


233 


106 


87 


14.1(1.7) 


1.7(0.5) 


G053.81-00.00.. 


MMl 


1.9 


74 


0.14 


6.9 


3.6 


28 


36 


28 


12.4(1.4) 


2.1(0.5) 



Notes. The first two columns give the HEC name and the millimeter clump number, the following columns represent (in order of appearance): 
kinematic distance, visual extinction based on A'h, from the 1 .2 mm emission, size at FWHM, hydrogen column density, hydrogen volume density, 
mass, mass within 0.25 pc diameter, virial mass, NH3 rotational temperature, and NH3 column density. 



The extinction mass is independent of temperature, but it still 
depends on the distance via the cloud surface A. We derive ex- 
tinction masses for the clouds from ~30 to 6500 M© (Table|2]l. 



5.2. Column densities and visual extinction 



As mm emission from cool clouds is usually optically thin, the 
column density and the mass of a cloud are well sampled by 



the observed flux density. The beam averaged column density is 
given bv Motte et al. (2007); 



1.2mm 



Qoth/^/Ci .2mmfi 1 .2mm(rdust) 



(7) 



where 2mm peak flux, Q the beam solid angle, a-] 2mm the 
dust opacity at 1.2 mm per unit mass column density, assuming 
a gas-to-dust ratio of 100, Bi.2mm(rdii.st) the (full) Planck func- 
tion at the dust temperature, and me and ji as defined before. 
The column density depends strongly on the dust properties: the 
dust opacity, k,\, and the emission coefficient, are related as 
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Fig. 6. Number distribution of 1.2 mm clumps with size deter- 
mined from Gaussian fits and converted to physical units using 
the near kinematic distance. 



HEC 0016,93 + 00,2-1 
MMl 
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LSR Velocity [km s '] 
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MM2- 
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Fig. 8. The NH3(1,1), NH3(2,2) and NH3(3,3) spectra taken 
with the Effelsberg 100m telescope toward two HEC clumps. 
The green lines indicate the theoretical intensity of the hyperfine 
lines given the intenstity of the main Une. 



*'o(;^) ^- We considered two different opacities; /fi.2mm 



Fig. 4. Examples of extinction maps (red dashed contour, start- 
ing at Av=21 mag and increasing by 10 magnitudes, beam 54") 

ontopof the 1.2 mm emission observed with MAMBO-2 (beam 0.4 cm after [H ildebrand (1983) and /ci. 2mm = 1.0 cmV^ 

10."5) in greyscale and solid black contours. taken from lOssenkopf & Hennina d 19940 . Table 1 column 6, 
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DV NH3(1,1) [km s ] 

Fig. 9. Number distribution of the NH3(1,1) line widths of the 
1.2 mm clumps. 



for dust grains with thin ice mantles. The emission coefficient 
was kept at P - 2, the advocated value for cold dust clumps 
et al.ll2006[). The dii st opacities differ by a factor 2.5, mean- 
ins the K of Hild ebrandl dl983h results in 2.5 times larger col - 
umn densities than the opacity of lOssenkop f & He nning (|1994|). 
In this work we used the dust opacitv of lOssenkop f & HenningI 
(Il994l) . A-i 2mm - l.Ocm^g"'. For the dust temperature we as- 
sumed the rotational temperature derived from ammonia (see 
Section I5.5l l. Since the clumps have high densities, collisions 
will dominate over radiative processes and the temperature ex- 
change between dust and gas will be efficient. In sources without 
NH3 detection we assumed a dust temperature of 16 K, which 
was the average NH3 rotational temperature. For clumps which 
have an embedded protostar the NH3 rotational temperature, de- 
rived on a 40" scale, might underestimate the dust temperature 
leading to an overestimation of the derived masses and column 
densities. 

We found clumps with column densities of the order of 
10^^ - 10^^ cm"^ (Table |7]). We derived corresponding peak vi- 
sual extinction values by applying Eq. |5] The peak fluxes, col- 
umn densities and peak visual extinction of Table |7] correspond 
to the positions listed in Table [3] The clumps have peaks in Ay 
from 31 to 280 mag, with an average value of 75 mag. These 
values are much higher than the mean visual extinction reached 
by the extinction method for this selected sample, which is be- 
tween 16 and 47 mag. This is expected because the peak val- 
ues are larger than the mean and with the limited resolution of 
54" higher extinction peaks (as found for example in IRDCs) 
are missed. Additionally, also the temperature might play a role, 
since the derivation of the Ay from the mm emission depends on 
temperature while the extinction method does not. 



5.3. Masses from 1.2 mm emission 

The clump mass can be derived fro m the 1 .2 mm emission by 
(iHildebrandl 1 983t iJMotte et"ani2007h : 



M 



1.2mm ~ 



F' 

1.2mm 



'<'l.2mm^l.2mm(7'dust) 



(8) 



where F[ j^^^ is the integrated flux density of the clump, d the 
distance, and Ki2mm, fii.2mm, and Tdust as defined before. The 
measured flux density can be contaminated by free-free emis- 
sion if an ionizing source is present. Since our sample contained 
very young objects, we can neglect such contamination. The 
derived masses depend strongly on temperature and distance; 



a 10% nearer distance decreases the mass by 20%, a 10% de- 
crease in temperature increases the mass by 17%. Nevertheless, 
the largest uncertainty in the mass derivation is caused by the 
dust properties, as described in the previous section. 

The majority of the clumps have masses between 12 and 
700 Mq and are located at distances between 2 and 5kpc (see 
the squares in the left panel of Fig. Q and Table |7]i. We find no 
low-mass clumps (M < 100 Mq) further than 4kpc, which is an 
bias of our extinction method (discussed in Sect. 16.41 ). Almost all 
high mass clumps (M > lOOM©) are located at distances larger 
than 4kpc. 

A second method was used as a comparison for the clump 
masses to check if the source finding algorithm was biased to a 
source size (and hence clump mass). Based on the (near) kine- 
matical distances, we defined circles of 0.25 pc diameter for each 
source. Then, we derived the integrated flux, Fo.25pc, and the 
mass, Mo.25pc, for the region within this circle (see Tables [3] Q 
and|4]i. The crosses in the left panel of Fig. Q show Mo.25pc as a 
function of the distance. The clump masses within 0.25 pc show 
a similar behavior as the clump masses determined by Gaussian 
fits. On a larger scale, the cloud mass was derived from the 
1 .2 mm emission according to Eq. [8] using the integrated flux 
down to 3cr. The cloud diameters are of order ~90", so after 
resampling the extinction map, it was possible to compare the 
cloud masses derived by extinction with the masses from the 
1 .2 mm emission. We find the extinction masses to be larger by 
a factor ~ 1 .3 than the masses derived from the dust continuum 
maps (see Fig.fTOb. This is expected, since the bolometer filters 
out large-scale structures by the sky noise subtraction and chop- 
ping. The cloud masses derived from the 1.2 mm emission and 
the extinction are listed in Table |2] 




10 10 1 

1.2 mm continuum mass [Mg 



Fig. 10. The cloud masses estimated from the color excess versus 
the masses estimated from the 1 .2 mm continuum maps. The red 
dashed line is a help line indicating the trend in which both mass 
estimates are equal. 



Finally, we estimat ed the volume-aver aged gas density of the 
clumps, riH,, following rMotte et al] (l2007h : 



M 



(9) 



where M is the clump mass, R the clump radius given by the 
geometrical mean of the semi major and semi minor axis from 
the Gaussian fit. The average gas density is ~ 2 x 10^ cm"^ and 
the individual results are given in Table |7] 



Rygl et al.: Initial phases of massive star formation in high infrared extinction clouds 



17 



o 



10' 

l«g(Mi,2„,n,) 



[Me] 



Fig. 11. Virial clump masses versus the 1.2 mm continuum 
clump masses. The red dashed line indicates the trend where the 
virial parameter is unity. 

5.4. Virial masses 

The thermal line widths are usually one order of magnitude 
smaller than the observed NH3 line widths, which indicates 
that the observed line widths are dominated by turbulence, 
hence contain information of the average kinetic energy within 
a clump. 

Given an optically thin line width and a c lump radius, the 
virial mass, can be calculated dMacLaren et al.l 

Il988h . Here cr is the three dimensional root-mean-square veloc- 
ity, R the clump radius and ki a density distribution constant. For 
a constant density distribution ki - 5/3. After the conversion of 
cr to the observable FW HM line width Av, cr^ = (3/81n2)Av^ 
dRohlfs & Wilsonll200l . the virial mass can be written as: 



/ Av x2/ R N 
M„„ = 210xf: .) (_) 



^kms 



[Mo]. 



The resultant virial masses. My 



are listed for each 
in Table [71 The virial parameter is defined as 



a 



(10) 



clump 

M 



(iBertoldi & McKeell992h . For a < 1 , the clumps are dominated 
by gravity, however with a » 1, the clumps are confined by 
surface pressure and self-gravity is unimportant. Most clumps in 
high extinction clouds have clump masses larger than their virial 
masses (see Fig. fTTT i. indicating that possibly most clumps are 
dominated by gravity and are collapsing. 

5.5. Temperatures 

The am monia molecule, described thoroughly bv lHo & TownesI 



1983h. is often used as a molecular cloud thermometer 



(iDanbv et al.| ll988). ts energy levels are parameterized by the 
/ and K quantum numbers, measuring, respectively, the mag- 
nitude of the total angular momentum and its component along 
the symmetry axis. Each set of rotational transitions is aiTanged 
into so-called /T-ladders, levels of fixed /T-value. From the sym- 
metry of the electric dipole moment of the molecule, all dipole 
transitions with nonzero AK are forbidden, meaning that the K- 
ladders are independent of each other The lowest transitions 
of each /T-ladder are metastable, and can be excited via colli- 
sions. Additionally, NH3 also undergoes vibrational motion from 
the tunneling of the nitrogen atom through the hydrogen plane, 
which splits the rotational energy levels into inversion doublets. 



The inversion transitions are further split into five 
quadrupole hyperfine lines, which allow the calculation of the 
optical depth (see the hyperfine lines in the NH3(1,1) transition 
in Fig. [8). Typical values of the optical depth were between 1 and 
3, showing that ammonia is optically thick in most cases. With 
the optical depth known, the rotational temperature follows from 
the ratio of the peak intensities of NH3(1,1) and NH3(2,2) lines 
after .Man gum et al. (1992): 



r r 0.283AV22 r 

* main^^ 1 1 



' 22 



(i-O 



[K],(ll) 



where Tmain is the NH3(1,1) main group optical depth, Avn, 
Av22, Til, T22 are the line widths and the peak intensities of the 
NH3(1,1) and NH3(2,2) lines, respectively. 

The line width ratio in this equation is debatable. 
Observations show that the NH3(1,1) and NH3(2,2) line widths 
are dominated by turbulence. If one assumes that both transi- 
tions arise from the same region, they experience the same tur- 
bulence, and the line widths are equal, AV22 - Avu. In this 
case, the linewidth ratio drops out of th e equa tion, which be- 
comes the equation from iHo & TownesI (Il983l) . We observed 
slightly larger NH3(2,2) line widths for several clumps, probably 
since this transition has a higher energy and is therefore more 
sensitive to warmer and more turbulent regions. The NH3(2,2) 
line widths are therefore more sensitive to the peak tempera- 
ture, determined by the turbulent outliers. We considered it rea- 
sonable to use Av22 =Avii, since it would give an 'average' of 
the NH3(2,2) turbulence, and hence an 'average' rotational tem- 
perature of that region. The tw o methods yield very similar re- 
sults; the lHo & TownesI ( Il983h temperatures wer e slightly lower 

Arrot - IK) than when using the formula of iMangum et al.l 

1992h . 

The clumps are on average cold ~16K. The temper- 
atures range from 10 to 25 K. For such low tempera- 
tures, the kinetic temperature is well approxi mated by T^ot 
dWalmslev & Ungerechtslll983l:lDanbv etal][T98 8). 

We derived colu mn densities of the NH3(1,1) line (after 
iMangum etaT]ll992l) : 



A^(i,i) = 6.60 X lO'^^AviiT— [cm-2 ] 



Vll 



(12) 



From t his we estimat ed the total ammonia column density (fol- 
lowing jLi etal]|2003h : 



(1,1) 



^ 1 .23.K 5 /-41.2x 
l + -exp(— )+-exp(— -) 



14 . 
+ y exp( 



Tvot 3 
-99.4 



Tro 



[cm-2]. 



(13) 



The averaged ammonia column density for the clumps is 3.1 x 



10 cm" . The NH3 rotational temperatures and column densi- 
ties for all clumps are given in Table |7] together with the other 
parameters derived in this Section. 



6. Discussion 

6.1. Higli extinction structures 

Thanks to the high resolution Spitzer IRAC data, the ex- 
tinction method allows to follow the mass distribution from 
Galactic size-scales down to single clouds. High extinction 
traces Galactic structure (see the Galactic distribution of the 
color e xcess in Fig. [2ll, similar to the CO survey of the Galactic 
plane dPame et alJ [l987i) and the dust continuum surveys. 
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ATLASGA L at 850 ([Schuller et al.ll2()()9t) and BOLOCAM 
at 1.1 mm ("Roso lowskv et alJ l2009h . In addition, the distribu- 
tion of 6.7 GHz Class II methanol masers traces high-mass star- 
formi ng regions and, thus. Galactic structure (iPestalozzi et al.l 
l2005h . All these surveys, except for the methanol masers, peak 
toward the Galactic Center The extinction maps miss the inner 
1° around the Galactic Center, since the Spitzer IRAC data were 
not publicly available for this region at the time. However, a ris- 
ing trend toward the Galactic Center was observed. The extinc- 
tion distribution indicates that the column densities are higher 
towards the fourth quadrant, 0° > / > -60°. We find no evi- 
dence for this from the CO distributi on, nor from th e methanol 
masers. Only the ATLASGAL survey (ISchuller et al.t200 9) hints 
at a similar distribution, but needs to be extended in longitude 
range for giving conclusive evidence. 

The next eye-catching features in the extinction distribution 
were two bumps, peaking roughly around longitudes of ~ 40° 
and -30° (see Fig. JUl, which were confi rmed by the methanol 
maser distribution (IPestalozzi et al.ll2005h . We identified several 
spiral arms by comparing our results to the Galactic models of 
[Vallee (2008). In Fig.|2]each spiral arm traced by a maximum in 
the extinction distribution is marked. We can trace the tangen- 
tial point of the Sagittarius-Carina arm at / ~ 50°, the tangen- 
tial points of the Scutum-Crux arm at / ~ 35° and / ~ -45°, 
the beginning of the Norma-Cygnus arm on the near end of the 
Galactic bar I ~ 10°, the beginning of the Sagittarius-Carina 
arm at the far end of the Galactic bar I ~ -5°, the beginning 
of the Perseus arm I ~ -20°, and the tangential point of the 
Norma-Cygnus arm at I ~ -30°. Several of these features are 
also seen in the ATLASGAL survey, and to a lesser extent in the 
CO distribution. Hence, the total extinction distribution seems 
to agree with results of previous studies, and tracks large scale 
mass structures. 

The combination of submillimeter dust continuum emission 
and extinction maps extends the size scales down to clump 
sizes of fractions of parsecs. In Table [8] the mean masses, radii, 
and volume densities for different size-scales of complexes, 
clouds and clumps are compared. Complexes are the low den- 
sity (Ay > 21 mag or A^h^ > 2 x 10^^ cm"^) regions surround- 
ing the cloud. In general, the masses and sizes decrease towards 
the smaller scales, while the volume density increases as ex- 
pected. The cloud masses derived from the extinction are higher 
than from the mm emission, as discussed in Sect. 15.31 The com- 
plexes and the clouds are confined by the pressure of the surroud- 
ning medium rather than by gravity. This is based on their virial 
masses, where we used for the line wi dth estimate the '-^CO data 
from the Galactic Ring Survey (GRS) d Jackson et al.ll2006l) . The 
'^CO line is a probe of the low density material and with the 
GRS resolution of 47" the line width is represenative for a cloud- 
scale average. Most of the clumps inside the clouds are, however, 
bound objects for which Mvir < Mi 2mm- 

6.2. Evolutionary sequence 

Our observations suggest different classes of clouds and in the 
following we discuss their possible relation to diff'erent evolu- 
tionary stages. Clouds, which either have no clumps or clumps 
of which the peak in the mm emission is less than twice the mean 
emission in the cloud, were defined as diff'use clouds - this def- 
inition is based on the cloud morphology and therefore difffer- 
ent f rom the "classical diffuse clouds" defined by Ay < 1 mag 
(ISno w & McCall 2006). Clouds which contain clumps with a 
higher contrast, above twice the mean cloud emission, were con- 
sidered to be peaked clouds. If the peak flux was above thrice 



Table 8. Mean properties of high extinction complexes, clouds 
and clumps 



Size scale 


Method 


Mass 


Radius 


Vol. den. 






(Mo) 


(pc) 


(cm-3) 


complex^ 


extinction 


4000 


1.4 


6x10' 


clouds 


extinction 


910 


0.7 


IxlO* 


clouds 


1.2 mm emission^ 


700 


0.7 


9x10^ 


clumps 


1.2 mm emission 


130 


0.15 


2x10' 



° The high extinction complexes should be treated as an indication 
since they have not so strictly defined boundaries as clouds and clumps. 
* smoothed maps by 20" Gaussian 
unsmoothed maps 



the mean emission and there were two or more clumps the cloud 
was classified as a multiply peaked cloud. The ratio of the clump 
peak emission to the mean emission of the cloud is given in Table 
|2] together with the classification. The mean physical properties 
such as temperature, masses, column densities and line widths 
are put together in Table|9]for the three classes. 

Diffuse clouds are the most likely candidates to form, or 
harbor, starless clumps, which are expected to be cold, more 
extended and more massive. The clump, or cloud, consists of 
low column density material and contains no massive compact 
object. Indeed, for very diffuse clouds, the source find algo- 
rithm yielded few clumps. On the smoothed maps, the algo- 
rithm returned values which were close to cloud size-scales with- 
out much substructure. This supports the idea that these diffuse 
clouds represent the earliest stage in which few condensations 
have formed and gravitational collapse has not yet started. More 
evidence for the young nature of diffuse clouds was found in the 
24 fim MIPS data: generally the (smoothed) dust emission fol- 
lowed the 24 jum-dark regions and only a few 24 yum sources of 
~100MJy sr"', located within 20" of mm peak, were found to- 
ward the clouds. In two cases 24 jum source > 100 MJy sr"' were 
found the edge of the cloud, at ~ 1 ' from the mm peak. Examples 
are G012.73-00.58, G013.28-00.34, and G034.85+00.43 (see 
Fig.©. 

As the first clumps become more compact and their mm peak 
flux rises above twice the mean emission of cloud, in the peaked 
cloud stage, the clump is accreting material from a lower density 
reservoir or envelope surrounding it, and consequently the tem- 
perature and turbulence rise. Since star formation is a process 
which is dynamic and includes feedback and triggering, akeady 
in this early stage a cloud can contain more than one clump, 
though possibly the other peaks barely stick out above the mean 
mm emission of the cloud. The peaked clouds are generally 
dark at 24 fim, except for a bright 24 fim source located within 
20" of the mm peak (see G013.91-00.51, G016.93-h00.24, and 
G053.8 1-00.00 in Fig.©. The peak flux of the 24 fim sources is 
on average ~ 250 Mjy sr ' , which is higher than what was found 
for the diffuse clouds. 

After more accretion from the reservoir, the clump will have 
a much higher mm peak flux compared to the initial stage, reach- 
ing at least trice the mean emission of the cloud. The clump mass 
will have decreased toward the more evolved stage, since not 
all of the mass reservoir will not be accreted onto the clump, 
as the star formation efficiency is apparently less then 100%. 
The column density, temperature and turbulence will increase 
even further than in the peaked cloud stage. While in the peaked 
stage, there were few clouds which had multiple clumps, in the 
more evolved stage one expects to find more clouds with mul- 
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tiple clumps. Most of these clumps should be well above the 
mean emission of the cloud. In these multiply peaked clouds, 
the clumps will be in different cloud stages of evolution depend- 
ing on the initial conditions of the clump. Also not all clumps 
will form high-mass stars, so they will have different properties. 
And indeed, several clumps in the multiply peaked clouds show 
bright, > 500 Mjy sr"', 24 jum emission, indicating a protostar, 
within 10" or less of the mm peak, while other clumps are in- 
frared dark. Such clouds with clumps in various phases are not 
completely infrared dark, since the dust is locally heated from 
the already hot HMPO or UCHii region. Examples for multiply 
peaked clumps are shown in Fig. |5) G014. 63-00. 57, with the 
very bright component MMl and the cloud G018.26-00.24 with 
five clumps. 

Additionally, one can check the evolutionary stage by 
searching for signs of star formation such as water masers. Class 
II methanol masers, and shocks. Table |9] also contains entries 
for such star formation tracers. Water masers are thought to be 
caused by outfl ows during (low and high-mass) star formation 
(lMentenlll996l) . The Class II methanol masers are uniquely as- 
sociated with high-mass star formation, and are usually found 
prior or coexistent with an H ii region (Menten 1991; Ellingse^ 
I2OO6 ; Pestalozzi et al. 2007). While the water masers are found 
towards both the peaked and the multiply peaked clouds, the 
methanol masers are only found toward the latter class. Another 
indication for star formation comes from gas excited by shocks, 
such as H2(v=0-0), 5(9,10,11) or CO(v=l-0) gas. The IRAC 
4.5 fim band contains both these transitions, and can therefore 
be used to search for ongoing star formation. The close asso- 
ciation of the ext ended 4.5 //m obje c ts w ith Class II methanol 
masers, found bv Cvg anowski etal] 1*2008), suggests that they 
might trace ongoing massive star formation. We find two multi- 
ple peaked clouds in our sample that contain such an extended 
4.5 fim source, of which one also contains Class II methanol 
masers. 

The peaked and multiply peaked clouds show active star for- 
mation, while the diffuse clouds do not. Possibly, these clouds 
might not form stars or will not have column densities high 
enough to form massive stars. However, some, for example 
0013. 28-00. 34, seem to be forming clumps which might evolve 
into active star-forming clumps. Table |9] shows that the clump 
temperatures are very similar, especially in the case of peaked 
and the multiply peaked clouds; the diffuse clouds stand out with 
their temperature of 13.5 + 1.5 K against the average tempera- 
ture of all clouds 16 + 1.8 K. As star formation produces stars 
in a large range of masses, so cluster-forming regions or clumps 
within one cloud are expected to have different column densities 
and temperatures. Given the range of column densities for the 
clumps 3 - 30 X 10^^ cm"^ not all clumps might evolve into clus- 
ters with massive stars, even if mass accretion is still continued 
from the low density material surrounding the clump. 



6.3. Comments on individual sources 

We performed a systematic search using the SIMBAD 
Astronomical Database of our cloud to find signs of star forma- 
tion such as maser emission, Hii regions and IRAS sources. We 
discuss here several diffuse, peaked and multiply peaked high 
extinction clouds, which are representative examples of typical 
clouds in our sample. Additionally, several unique and interest- 
ing clouds are mentioned. 



Table 9. Mean physical properties of the clumps and star for- 
mation indicators in diffuse, peaked and multiply peaked clouds. 





Diffuse 


Peaked 


Multiply peaked 


Mass (Mo) 


230^ 


185 


150 


Si 7P ( np^ 


0.5~ 


0.36 


0.30 


TVnh, (102-cm-2) 


4.5" 


5.5 


8.5 


r,-o, (K) 


13.5 


1 1 
1 J. / 


n s 
i / .J 


Av (kms"') 


1.2 


1.4 


1.6 


IRDC corr.^ (within 2') 


83% 


50% 


70% 


24/jm^ (within 1') 


17% 


83% 


100% 


Ext. 4.5pm^ (within 1') 


none 


none 


29% 


Masers_^ (within 1') 


none 


H2O 25% 


H2O 57% 






CH3OH 29% 


" As most diffuse clouds contained 


no clumps, 


we derived the sizes 



from the smoothed bolometer maps of the diffuse clouds. These are 
therefore not properties of the clumps in diffuse clouds, but properties 
of the diffuse clouds themselves. 
* Catalogs of lSimon"et 

2] ( I2006j) and lSimon et"^ ( l2006bh . 

> lOOMJy sr-' 

Ex tended 4.5/jm sources from the catalog of ICveanowski et al.l 
J200a) 

' H 2O: iJaffe et al.l jl98lh and this work, CH^OH: ISzvmczak et"an 
l l2OO0l) 



6.3.1. G01 3.28-00.34 

The cloud GO 13. 28-00. 34, with its one centrally located clump 
MMl, is a good example of a cloud in a very early phase of 
protocluster formation. This clump has a low column density 
3 X 10^^ cm"^ and contains less than 5% of the total cloud mass. 
Likely, this clump did not yet start to collapse, since from the 
virial mass we know that the clump is confined by gas pressure. 
Additionally, the clump contains no 24 /vm emission, which is a 
further evidence for stage prior to star formation. The moderate 
NH3 line widths of 1 .6 km s ' indicate a source of turbulence or 
large scale motions, however, without an NH3 map of the cloud 
it is impossible to isolate the source for the turbulence. 

6.3.2. G013.91-00.51 

In the cloud GOO. 91-00. 51 the 1.2 mm emission follows an in- 
frared dark patch. The one clump in this cloud which is three 
times above the mean emission is located nearby two infrared 
emission peaks and can be possibly triggered by the star forma- 
tion in these infrared bright sites. The clump is very massive, 
145 Mq, and cold. The virial mass is 50 M© suggesting a col- 
lapse. 

6.3.3. G01 4.63-00.57 

Cloud G014. 63-00. 57 is nested around two bright infrared 
sources. The 1 .2 mm peak fluxes are high, especially for clumps 
nearby the infrared sources. The brightest clump, MMl, con- 
tains a compact bright infrared source, I RAS 18164-1631, with 
a flux density at 25yum Fii^,^ = 23 Jy (iHelou & Walkel^ll988h 
and has wide NH3 line widths of 1 .8 km s~' indicating large tur- 
bulenc e. Additionally, w ater masers were detected tow ards this 
clump (Jaff^e etal." 198 1") and 'Cvganowski et al.1 (12008') showed 
that there is extended 4.5 jim emission - a further indication for 
outflows. All evidence points to a state of ongoing star formation 
for the MMl clump. The other clumps of G014.63-00.57 are 
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R.A. (jaooo) 

Fig. 12. A Spitzer GLIMPSE 8.0 yiun image of cloud 
G017. 19+00.81 overlaid with the 1.2 mm emission in white 
contours starting at 2cr (cr - 0.015 Jy beam"') and increas- 
ing by V2. All clumps are at the same kinematic distance of 
2.3 + 0.2 kpc. The dashed circles guide the eye to the arc-like 
structure, which is infrared dark in the north-east and infrared 
bright in the south west. 



still infrared dark. The MMl and MM2 clumps have the theoret- 
ically required column density for massive star formation (see 
Sect.lEHl. 

6.3.4. G01 6.93+00.24 

In G016. 93+00. 24 the 1.2 mm emission follows the infrared 
dark filament. In the west end of the filament, the 1 .2 mm peak 
MMl, is a few arc seconds offset from a infrared bright emission 
peak. The NH3 line widths are relatively narrow, 0.9 km s and 
the temperature is low, 14 K, which suggests that the clump is in 
a very early phase with little turbulence. 

6.3.5. G017.19+00.81 

The multiple peaked cloud G017. 19+00.81 harbors several 
clumps, which are all at a kinematic distance of 2.3 ± 0.2 kpc, 
arranged in an infrared dark arc (Figs. ISl andfT2l). Such a mor- 
phology might point to triggered star formation. The location of 
this cloud, on the edge of the bright Hii region Ml 6, the Eagle 
nebula, suggests that we are observing new massive star-forming 
clumps in an early stage of evolution, whose creation might 
be linked to the already ongoing star formation in this region. 
G017. 19+00.81 differs from the classical triggered star forma- 
tion scenario described by Zavagno et al. (2005). First, the ab- 
sence of an Hii region in the center of the arc; no radio emission 
was d etected either in the 20 cm MAGPIS survey (' Helfand et al.l 
I2OO6I) . the 21cm V LA survey (Condon etal. 19981), or in the 
6 cm Parkes survey (iHavnes et al. 1979) - the object was not 
covered by the 5 cm VLA survey (Becker et al. 1994). Second, 
an arc triggered by an Hn region is expected to be infrared bright, 
whereas G017. 19+00.8 1's arc is partially infrared dark. 



The clumps in the cloud have similar temperatures and sizes. 
The upper two left clumps, MM2 and MM3, represent most 
of the mass in this cloud, whereas the southern MM4, shows 
the largest amount of turbulence. Toward the MM2 clump we 
detected water masers, within a beam of 40". It is coincident 
with a bright 24 jum source indicating a protostar. MM2 seems 
to be an evolved clump, where star formation is just begin- 
ning. Additionally, the column density of this clump theoreti- 
cally allows high-mass star formation. The MM4 clump is lo- 
cated nearby an Hii region, which might explain the wide line 
width of 2.8 km s"' and high degree of turbulence. Nearby (1') 
there is an infrared-bright source, IRAS 18164-1340, which 
is only detected at long mid-infrared w avelengths of lOOjum 
ioo;,m = 770 Jv. iHelou & WalkedlT988l) We performed a more 
detailed study of this region; a molecular line survey with IRAM 
30m and an interferometer study with the SMA (Rygl et al. , in 
preparation). 

6.3.6. G01 8.26-00.24 

GOl 8.26-00.24 is a cloud with five 1.2 mm peaks located north- 
ward of an infrared bright region with several Hn regions 
and water maser detections. There are several infrared sources 
present in the cloud, however, the most mm emission peaks 
are infrared dark. MM3 is an exception, here the infrared peak 
is just a few arc seconds from the mm peak. MMl has the 
largest separation from the infrared bright region. This clump 
has the highest column density and the narrowest line width 
of the five clumps in this cloud. The overall cloud shape sug- 
gests that the gas and dust was collected by a driving force 
south-west of the cloud, possibly the infrared bright source at 
(a = 18''25™02.''8,5 = -13°09'30", J2000). 

6.3.7. G022.06+00.21 

G022. 06+00. 21 contains one clump with a bright mm 
peak, MMl, and second clump with a weaker peak, 
MM2. MMl is coincident with IRAS 18278-0936 
(f25;.m = 12 Jy, 'Helou & Walkei'119881), has 6.7 GHz methanol 
masers (ISzvm czak et al. 2000), extended 4.5 //m emission 
(Cvganowski et al.l 120081) . and has a sufficient column density 
for high-mass star formation. Additionally, we detected water 
masers towards MMl, within a beam of 40". MMl is a rela- 
tively warm clump, 25 K, with ongoing star formation suggested 
by the masers and outflow tracers. MM2 is more quiescent and 
coincident with a weak infrared source. 

6.3.8. G024.37-00.15 

G024. 37-00. 15 contains two clumps, one infrared bright, MMl, 
and one infrared dark, MM2. MMl coincide s with the bright 
IRAS source 18337-0743 {F^^^^ = 47 Jy, iHelou & Walked 
1988), and has a wide linewidth 2.1 kms"'. The MM2 core 
is in an earlier state of star formation, exhibiting SiO emission 
dBeuther & Sridharanll2007h and water masers. 

6.3.9. G024.61-00.33 

In G024. 6 1-00.33, the strongest mm emission peak, MMl, 
is coincident with IRAS 18346-0734 (Fas^,™ = 25 Jy, 
Helou & Walke r 1988) and nearby, separated by ~ 30", 
ISzvmczak et al.l ( l2000l) detected 6.7 GHz methanol masers. The 
second clump, MM2, is infrared dark. 
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6.3.10. G024.94-00.15 

The mm emission follows closely the infrared dark filament in 
cloud G024.94-00. 15. The clumps, MMl and MM2, are located 
near weak infrared sources. We found water maser emission to- 
ward MMl, within a beam of 40", indicating early phases of star 
formation. 



Given the size of our bolometer map, we identified the same 
clumps within this cloud. Taking into account different analysis 
metho ds, we find similar fluxes and masses as .Rathborne et all 
(l2006h . 



6.3.1 1 . G030.90+00.00 A, B, C, and D 

G030. 90+00. 00 consists of several high extinction clouds at very 
different distances (Fig.|5j. It is located at the tangential point of 
the Scutum-Crux spiral arm (Vlsr ~ 105 km s"'), where it over- 
lags with the Sagittarius-Carina (Vlsr ~ 35 and ~ 75kms"') 
( IVallee.,2008i) . The Scutum-Crux arm connects at this longitude 
of 30° with the Galactic bar, a region with highl y shocked gas 
and explosive star formation (iGarzon et al. featuring ad- 

ditionally the mini-starburst region W43. We found one clump, 
MM2, which, considering the kinematic distance at 7.2 kpc, is 
located on the end of the bar. 

6.3.12. G034.77-00.81 

G034.77-00.81 is a cloud with very diff'use 1.2 mm emission lo- 
cated in a infrared dark part of the sky. By smoothing with a 
20" Gaussian the extended emission became more evident. We 
detected weak NH3(1,1) emission toward the center of the ex- 
tended emission. The ammonia emission indicates that there is a 
dense region in this cloud, however, it has too low column den- 
sity or is too extended to be detected by MAMBO-2. Such clouds 
as G034.77-00.81 could be the most early stages of within our 
sample. 

6.3.13. G035.49-00.30 A and B 

Cloud G035. 49-00. 30 A contains one clump, MMl, coincident 
with a bright infrared source. In the nearby G035.49-00.30 B we 
found five clumps located in an elongated infrared dark cloud 
with several weak infrared sources. We detected water maser 
emission towards MMl and MM3, within a beam of 40". MMl 
is the most evolved clump, with ongoing star formation. MM3 
should be in a earlier stage of star formation, since the peak of 
the mm emission is infrared dark. 
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Fig. 13. The NH3(1,1) line width as a function of the NH3 rota- 
tional temperature. Linear regression finds a weak trend, coiTe- 
lation coefficient of 0.46, of increasing line width with tempera- 
ture. 
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Fig. 14. The NH3(1,1) line width as a function of the NH3 col- 
umn density. In general, larger line widths are found for higher 
NH3 column densities with a correlation coefficient of 0.53. 



6.4. Comparison with infrared dark clouds 

Recent studies have shown that IRDCs generally ha ve a filamen- 
tary shape and a high column density of ~ 10^^ cm"^ dCarev et al.l 
[T998„ 2000: Rathborne et al.. 200 6). IRDCs, which are very com- 
pact or have large aspect ratios, might not be detected in the ex- 
tinction maps because of the limited resolution. The extincion 
method is more sensitive to the more extended and therefore 
on average lower column density material than generally found 
in IRDCs. However, the follow-up study of the 1.2 mm emis- 
sion maps with MAMBO is sensitive enough to reach the typi- 
cal IRDC column densities, allowing a comparison between the 
properties of the high extinction clouds and the IRDCs. 

We find that almost 70% of the high extinction clouds have 
also IRDC characteristics using the Catalog of IRDC candi- 
d ates from [Si mon et al. (2006a) and the Spitzer IRDC catalog 
of iPeretto & Fullen (i2009). This overlap caused that we inciden- 
tally studied the same source as Rathborne et al. (2006), namely 
G035. 4 9-00.30 (listed as MS XDC G035.39-00.33 in the pa- 
per by IRathborne et al.l l2006l) . allowing a direct comparison. 



In general, the temperature and the level of turbulence (line 
withs) of the clumps in hi gh extinction clou ds was similar to 
values found for IRDCs jPiUai et all |2006|) . Figures [13] and 
[T4l show that the turbulence tends to increase with ammonia 
column density a nd temperature, which is similar for IRDCs 
(iPillai etalJl2"006h . The masses and volume averaged densities 
of the clumps in our sample are similar to the averaged values 
for clumps in IRDCs (referred to as "cores" in Rathborne et al] 
2006). However, at large distances the high extinction method is 
biased tow ard more massive clump s. Comparing with the IRDCs 
sample of IRathborne et all (l2006h our sample is missing low- 
mass clumps (Mciump < 100 Mq) at large distances (d > 4 kpc), 
see the left panel of Fig.|7] However, our goal, to search for high- 
mass molecular clouds, is not affected. 

Although, the majority of the high extinction clouds is in- 
frared dark, there are a number of infrared bright sources. 
Especially amongst the multiply peaked clouds, we find several 
cases of already ongoing star formation via detections of Hii re- 
gions, water and Class II methanol masers, and extended 4.5 //m 
emission. Since we did not require the high extinction clouds 



22 



Rygl et al.: Initial phases of massive star formation in high infrared extinction clouds 



to be infrared dark, the method is also sensitive to the evolved 
stages of star formation. In addition, the color excess selection 
is sensitive to embedded young stars with a color excess of their 
own. 

6.5. Comparison with more evolved objects 

Later stages of massive star formation are the HMPOs and the 
UCHii regions. Since these objects are usually still embedded in 
their natal cloud, they emit most strongly in the millimeter to 
infrared. Additionally, UCHii regions emit free-free emission, 
which is detectable at radio frequencies. Most red and bright 
IRAS sources are such embedded stages of evolved star (or clus- 
ter) formation (Wood & Churchwell 1989). Ammonia observa- 
tions of HMPOs and UCHii regions show that they are usu- 
ally w arm (T > 22 K, ISridharan et 10120021 IChurchwell et all 
1 19901) and very turbulent (Une widths of 2-3 km s"^). Water and 
Class II methanol masers are observed towards these regions 
(IWalsh e t al. 1998), indicating outflows and infrared radiation, 
respectively. 

The high extinction clumps have lower temperatures and 
narrower line widths. Many of the more evolved clumps in 
high extinction clouds have their peak of mm emission nearby 
bright infrared sources. In several cases the two emission peaks 
even coincide. For three such cases, G014. 63-00. 57 MMl, 
G017.19-h00.81 MM2, and G022.06-i-00.21 MMl, we find 
slightly higher rotational temperatures, respectively 18, 19 and 
25 K, and detections of water masers reminiscent of outflows in 
the early phase of star formation. 

The clump masses in the high extinction clouds are al- 
most one order magnitude lower than the masses of HMPOs 
tBeuther et a l. 2002). Likely, by removing the median emission 
before performing the source extraction, the low column density 
mass reservoir of the clump is subtracted resulting in less mas- 
sive clumps than in the HMPO studies. Also the HMPOs are on 
average larger with an average FWHM of 0.5 pc, than the sizes 
of our clumps (~0.3 pc). 

In most of our clumps the ammonia column densities were of 
an order of magnitude higher than those for HMPOs. As a clump 
forms protostellar objects, the temperature and line widths in- 
crease, while the clump will start losing its envelope, which 
will decrease the ammonia column density on clump size-scales. 
Similar behavior of temperature, line wi dth and ammoni a col- 
umn density is also observed for IRDCs dPillai et al.ll2006l) . 

6.6. Comparison witli tlieoretical models 

There are two competing models to exp lain the for- 
mation of massive stars, core accretion (iMcKee & TanI 
[2003'; 'Krumholz et al.' '2005', '2009^ and competitive accretion 
(Klessen et al. 1998; Bonnell et al. 2001; Bonnell & Bate 2002, 
I2OO6 ). They differ primarily in how the mass is collected which 
will ultimately make the massive star In the first model, core ac- 
cretion, a massive star forms in a massive core, hence the mass 
reservoir of the core determines the mass of the star. In the sec- 
ond theory, a clump of gas will fragment into many low-mass 
objects containing masses around the thermal Jeans mass. The 
accretion will continue in this N-body system, where the overall 
system potential funnels gas down to the center of the potential, 
to be accreted by the massive stars forming there. A powerful 
way to distinguish between these scenarios are better determina- 
tions of the massive cores within the clumps. For the core accre- 
tion model these cores should be distributed similar to a stellar 



mass function. The degree of mass segregation of the cores in 
a clump and the limit of fragmentation will test the competitive 
accretion model. Such massive core studies require a resolution 
much higher than this study; they have to resolve the clumps 
into cores. Several studies on fragmentation were recently car- 
ried out using the Plateau du Bure or the SMA in t erferometers 
(Beuther et al.' 2007'; 'Rathb orne et"ani2007l 12008'; 'Zhan g et al.1 
2009; Swift 2009). They show that the 1.2 mm clumps fragment 
into smaller cores, but cannot tell at what level the fragmentation 
halts. With the advent of ALMA, the limit of fragmentation will 
finally come within observational reach. 

Simulations o f core accretion models by 

iKrumholz & McKe9 (l2008|) put forward a minium column 
density of 1 g cm"^ for massive star formation. For reference, 
this corresponds to 4750 Mopc'^ or - 2.7 x 10^^ cm"^. 
This limit would imply that only one clump from the high 
extinction clouds can form massive stars, namely G014.63- 
00.57 MMl. We note here that Ki'umholz & McKee (.2008) use 
a A:i.2mm ~ 0.4cm^g"', while we used /ci.2mm ~ l.Ocm^g"'. 
To compare this limit with our measured column densities we 
would have to scale our results by a factor 2.5, and in this case 
four clumps would theoretically be able to form massive stars. 
Of these four clumps only one, G014.63-00.57 MM2, is infrared 
dark. The other three clumps are infrared bright, indicating 
that they already commenced forming stars. It suggests that the 
theoretical column density limit is too high for the very initial 
stages of star cluster forming clumps. 

7. Summary 

We made infrared extinction maps of the first and fourth Galactic 
quadrant, compiling a catalog of compact high extinction fea- 
tures. We studied 25 high extinction clouds in detail, using 
IRAM 30m bolometer observations and pointed ammonia ob- 
servations with the Effelsberg Telescope. The main results of the 
paper can be summarized as: 

1. Using the average color excess of the Spitzer 3.6 - 4.5 yum 
IRAC bands, we are able to trace high density structures 
across the Galactic plane. The extinction selected clouds are 
found on distances between 1 and 7kpc. Most of them are 
concentrated in the range of 1 - 4kpc, and trace the same 
regions of the Galactic plane as IRDCs. 

2. The extinction method is more sensitive to large scale lower 
column density clouds than the bolometer observations of 
1 .2 mm emission, which resulted in a few non detections of 
clouds in the mm emission despite a mean visual extinction 
of ~30 magnitudes from the extinction map. In general, there 
is a good correlation between the cloud masses derived from 
the extinction maps and from the 1 .2 mm emission. 

3. From the 1.2 mm emission, we have found clumps with 
column densities of 3 - 30 x 10^^ cm"^ and masses of 
12 - 400 Mq. Evidently, not all clumps found by high ex- 
tinction will be able to form massive stars. It is expected for 
cluster-forming regions to produce low and high-mass stars, 
thus to have a wide range in masses. 

4. High extinction clouds contain a wide range of evolution- 
ary stages. 70% of the high extinction clouds are associated 
with infrared dark clouds, however several clouds show more 
evolved stages of ongoing (massive) star formation. Three 
different classes of clouds are proposed: 

- Diffuse clouds have no signs of masers, Hii regions or 
extended 4.5 //m objects, which are general signposts 
of star formation. The 1.2 mm emission shows no mm 
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peaks above twice the mean emission of the cloud. These 
clouds are cold {T ~ 13.5 K) and Une widths far above 
the thermal value mark the importance of turbulence. 
The observations suggest that the diffuse clouds are in 
an early phase, where the accretion of gas and dust into 
clumps might be on the verge of beginning or might 
never reach the necessary column density to form stars. 

- When the clumps manage to accrete more matter so that 
the mm emission exceeds twice the mean cloud emission 
they enter the next (peaked) stage. During the peaked 
stage already some clouds show more than one clump, 
which indicates that star formation can start in different 
regions in the same cloud at different times. 

- Multiply peaked clouds show many signs of ongoing 
star formation, such as detections of masers, H n re- 
gions or extended 4.5 jum objects. Several clumps are 
infrared bright and show generally slightly higher tem- 
peratures and turbulence than diffuse or peaked clouds. 
Clumps within a cloud are not necessary in the same 
state, and will hkely not all form massive stars. We find 
four clumps that satisfy the theoretical column density 
requirement for high-mass star formation. 
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Table 1. Catalog of high extinction clouds (HECs) with the peak color excess, £ < (3.6 jum - 4.5 jum) > 
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-07:40 


56 


0.32 


G333.31-00.36 


16:21:32.2 


-50:24 


41 


0.51 


G024.43-00.24 


18:36:42.2 


-07:40 


06 


0.31 


G333.47-00.15 


16:21:18.9 


-50:08 


51 


0.44 


G024.50+00.09 


18:35:39.2 


-07:27 


13 


0.27 


G333.49-00.24 


16:21:46.8 


-50:12 


19 


0.57 
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Table 1. continued. 



HEC name 


R.A. 


Dec. 




Color Excess 


HEC name 


R.A. 




Dec. 




Color Excess 






(h:m:s) 


^o. /. ft 


) 


(mag) 






(h:m:s) 






) 


(mag) 






(J2000) 


(J2000) 










(J2000) 




(J2000) 










1st quadrant clouds 












4th quadrant clouds 








G024 


61-00.33 


10. J / .ZL.O 


—K) 1 .jj 


08 





31 


G333 


60-00.22 


iD.z2.Uy 


Q 

y 


— jU.UO 


32 





44 


G024 


64+00.15 


1 Q.'2<./1 1 A 
io.jJ.4i.4 


on • 1 


22 





38 


G333 


66+00.37 


io.iy.ju 


Q 


/I Q.O 


47 





68 


G024 


82-00.11 


io.JU.J /.D 


— U / . i 


02 





26 


G333 


75-00.33 


10. Zj. i / 





— JU.U4 


46 





56 


G024 


94-00.15 


16. J /. iy.4 


— U/. lU 


34 





35 


G333 


76+00.35 


10. ZU. ZD 




— 4y.jj 


26 





54 


G025 


15-00.28 


lo.jo.uy.o 


— U/,Uz 


52 





32 


G334 


20-00.20 


10.Z4.4Z 


Z 


— 4y.jy 


56 





43 


G025 


63-00.12 


io.jo.jU.i 


— U0.-3Z 


55 





32 


G334 


45-00.24 


10. Zj. JO 


1 
i 


— 4y. J i 


03 





50 


G025 


79+00.81 


iO.JJ.ZO.O 


— Uj. JO 


57 





19 


G335 


06-00.42 


lo.zy.zz 


c 

J 


/I Q- 1 1 

— 4y. i i 


57 





77 


G028 


53+00.21 


1 Q ./I .on 1 


r\'2-/i Q 


20 





30 


G335 


25-00.30 


lo.zy. JO 


Z 


A C-<Q 

— 4o.jy 


04 





68 


G030 


48-00.38 


lo.4o.ZU.l 


— UZ.Zl 


16 





29 


G335 


28-00.13 


lo.zy.uu 


Q 


— 4o.jU 


43 





69 


G030 


62+00.18 


lo.40. J4.4 


— Ul.Jo 


31 





31 


G335 


44-00.23 


lO.jU.U/ 


U 


— 4o.4/ 


56 





82 


G030 


73+00.12 


io.4/.UU.U 


— Ui.j4 


22 





32 


G337 


15-00.39 


10. J /.40 


A 

4 


— 4/.3o 


49 





79 


G030 


90+00.00 


lo.4/.4j.l 


— Ul,4o 


24 





30 


G337 


45-00.40 


10. jy.uu 


Z 


—4 / .ZO 


07 





45 


G031 


05+00.27 


10.4/.UZ./ 


— Ul.ao 


08 





30 


G337 


50-00.19 


10.J6. 16 


n 

y 


— 4/.1J 


09 





48 


G034 


03-00.33 


16.J4. jo.y 




43 





30 


G337 


77-00.34 


10.4U.UU 


1 


— 4/.U0 


50 





53 


G034 


11+00.06 


lo.Jj.ZZ.U 


+U1 .U4 


37 





30 


G337 


93-00.51 


10.41 .Zj 


D 


/IT-HQ 
— 4 / .Uo 


51 





64 


G034 


34-00.90 


io.j /. iz.o 


+UU.JU 


40 





30 


G339 


26-00.41 


10.40.U1 


Q 


/I A-H/l 

— 4D.U4 


40 





53 


G034 


35-00.72 


lO.JO. JJ. / 


+UU.JO 


11 





34 


G339 


58-00.13 


10.40. UU 


1 
1 


— 4J.JO 


54 





52 


G034 


71-00.63 


lO.JO.JJ.O 


_i_ni • 1 T 


44 





41 


G339 


62-00.12 


10.40. Uj 


1 
1 


— 4j.j / 


07 





41 


G034 


77-00.81 


10. J /.4U.U 


1 ni . 1 c 


57 





39 


G340 


06-00.24 


1 /;./to. 1 A 
10.40. 14 


Q 

y 


/I ^-O 1 


32 





47 


G034 


85+00.43 


LO.JJ.Zi.O 


1 ni .C/1 
+Ui.j4 


07 





34 


G340 


26-00.24 


io.4o.jy 





iZ 


03 





57 


G034 


98+00.30 


1 Q.C/i .rv/r o 
i0.j4.U0.J 


+Ui.J / 


43 





31 


G340 


77-00.12 


lO.jU. 10 


K) 


A A -A A 

—44.44 


05 





44 


G035 


19-00.75 


1 Q.CQ. 1 'J Q 

lo.jo. 1 j.y 


+Ui.4U 


07 





59 


G340 


93-00.23 


lO.Jl.Zl 


A 

4 


AA-ACi 

— 44.4U 


36 





61 


G035 


49-00.30 


1 Q.<'7. 1 n yl 
16. J /. 1U.4 


+UZ.Uo 


27 





38 


G341 


12-00.42 


lO.jZ.jl 


Q 

y 


— 44.jy 


26 





60 


G036 


42-00.15 


10.J0./U.4 


+UJ.UZ 


04 





27 


G341 


12-00.42 


lO.jZ. J 1 





— 44. jy 


17 





61 


G037 


26+00.09 


io.jy.ui.y 




40 





36 


G341 


21-00.24 


lO.jZ.Zj 


J 


/I /I -OQ 


11 





47 


G037 


44+00.14 


io.jy.uy.o 


_i_f\A -CkA 

+U4.U4 


18 





28 


G342 


57+00.18 


lO.jj.Zl 


c 

J 


/I '2-r\Q 
—^■j.Kjy 


13 





56 


G037 


48+00.07 


i6.jy.zy.j 


_i_r\A .n^ 


03 





29 


G343 


40-00.33 


1 /.UU.Zl 


c 

J 


AO- Ad 

— 4Z.4y 


13 





52 


G037 


54+00.20 


16.jy.U6. / 


_i_(\A . 1 1 

+U4, 1 1 


44 





29 


G343 


75-00.16 


1 / .UU.46 


1 
1 


— 4Z.ZO 


19 





56 


G037 


65+00.12 


1 Q.CQ.OQ C 

lo.jy. jo.j 


_i_f\A • 1 ^ 


22 





29 


G343 


76-00.15 


1 / .uu.4y 


1 
i 


/I -0 < 

— 4Z.ZJ 


24 





56 


G038 


93-00.36 


iy.uj.4u.y 




19 





48 


G343 


84-00.08 


1 / .UU.4j 


1 
i 


/I • 1 

— 4Z . i y 


04 





66 


G044 


30+00.03 


1 Q. 1 0. 1 T f\ 
ly.LZ.l / .\J 


+ iU.U/ 


02 





31 


G344 


10-00.65 


1 / .U4.UJ 


y 


^2:27 


36 





81 


G046 


33-00.24 


ly.i /.uj.y 


+ 1 1.4 / 


23 





32 


G344 


21-00.61 


1 n.f\A • 1 A 
1 / .U4. 14 


Q 


-A2:2\ 


09 





44 


G048 


90-00.27 


iy.zz.u6.y 


+ 14. UZ 


57 





41 


G344 


99-00.23 


1 / .Uj. lU 


A 
4 
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38 





87 


G049 


39-00.31 
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,1/1 .0'7 

+ i4.Z / 


17 





34 


G345 


00-00.23 


1 /.Uj. IZ 


c 

J 


-41:29 


35 





50 


G049 


48-00.38 


iy./j.4u.o 


.1/1 .an 
+ 14. 


25 





60 


G345 


04-00.21 


1 /.Uj. Ij 


Z 


-41:26 


59 


1 


20 


G050 


06+00.06 


19:23:13.0 


+15:13 


28 





36 


G345 


26-00.04 


1 / .Uj. 14 


c 

J 


^1:10 


16 





80 


G050 


39-00.41 


19:25:34.4 


+ 15:17 


32 





35 


G345 


49+00.31 


1 / .U4.zy 


1 

1 


-40:46 


27 





64 


G053 


14+00.07 


19:29:18.0 


+ 17:56 
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46 


G345 


50+00.34 
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1 / .U4.Z4 
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J 


-40:44 


38 


1 


21 


G053 


21-00.09 


19:30:01.7 


+ 17:55 


28 





34 


G345 


67+00.34 


1 / .U4. JO 
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J 


-40:36 


44 





82 


G053 


24+00.06 


19:29:31.8 


+ 18:01 


14 





57 


G348 


18+00.47 


1 1 0. 1 n 
1 1 AZ. lU 
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4 


-38:31 


19 





95 


G053 


57+00.06 


19:30:12.6 


+18:18 


35 





37 


G350 


02-00.51 
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J 


-37:35 
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95 


G053 
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19:30:26.3 


+ 18:20 
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33 


G350 


52-00.36 
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J 
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79 


G053 


81-00.00 


19:30:54.4 


+ 18:29 
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41 


G350 


69-00.48 
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1 /.Zj.ZO 
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y 


-37:01 
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52 


G058 


48+00.42 


19:39:01.4 


+22:46 


46 





30 


G350 


94+00.75 
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/ 


-36:07 
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45 


G059 


63-00.18 


19:43:47.6 


+23:28 


45 
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94+00.66 
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z 


-36:09 


52 





83 


G059 


79+00.06 


19:43:13.2 


+23:44 
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29 


G350 


96+00.55 


17:19:58 


7 


-36:12 


56 
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4th quadrant clouds 








G351 


16+00.71 


17:19:53 


7 


-35:57 


29 





62 


G300 


91+00.88 


12:34:13.0 


-61:55 


39 


1 


00 


G351 


25+00.66 


17:20:23 


3 


-35:54 


39 





82 


G305 


36+00.19 


13:12:33.8 


-62:35 


12 





58 


G351 


44+00.66 


17:20:54 


7 


-35:45 


11 





92 


G309 


13-00.14 


13:45:14.1 


-62:21 


36 





52 


G351 


47-00.45 


17:25:30 


5 


-36:21 


33 





75 


G309 


37-00.12 


13:47:18.1 


-62:17 


35 





22 


G351 


52-00.56 


17:26:07 


7 


-36:22 


40 





40 


G309 


42-00.62 


13:48:37.0 


-62:46 


00 





51 


G351 


53-00.56 


17:26:06 


6 


-36:22 


05 


1 


30 


G310 


22+00.39 


13:53:21.6 


-61:36 


15 





66 


G351 


59-00.36 


17:25:26 


9 


-36:12 


41 





65 


G311 


57+00.31 


14:04:24.5 


-61:19 


45 





53 


G351 


78-00.54 


17:26:44 


9 


-36:09 


17 





58 


G311 


60+00.41 


14:04:25.5 


-61:13 


36 





50 


G351 


81+00.65 


17:21:59 


2 


-35:27 


34 





94 


G311 


95+00.15 


14:07:48.4 


-61:22 


56 





49 


G351 


96-00.27 


17:26:08 


4 


-35:51 


15 





71 


G312 


11+00.27 


14:08:47.1 


-61:13 


09 





47 
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Table 5. Observed ammonia parameters with uncertainties (in parenthesis) from hyperfine and Gaussian fits with CLASS 



HEC name 






NH3( 


1,1) 




NH3(2,2) 


NH3(3,3) 






Vlsr 




Av 


Tmain 


Tub 


Av 




Av 






(kms-') 


(K) 


(kms-') 




(K) 


(kms-') 


(K) 


(kms"') 


G012.73-00.58.. 


MMl 


6.48(0.01) 


1.3(0.3) 


0.7(0.1) 


4.1(0.4) 


0.2(0.1) 


0.8(0.1) 








MM2 


6.19(0.01) 


1.8(0.3) 


1.0(0.1) 


1.3(0.1) 


0.3(0.1) 


1.6(0.2) 






G013.28-00.34.. 


MMl 


41.30(0.06) 


3.2(0.7) 


1.6(0.2) 


1.8(0.5) 


1.3(0.3) 


1.6(0.2) 






G013.91-00.51.. 


MMl 


22.94(0.02) 


2.5(0.3) 


1.3(0.1) 


1.8(0.3) 










G013.97-00.45.. 


MMl 


19.76(0.02) 


1.4(0.2) 


2.3(0.1) 


0.7(0.1) 


0.5(0.1) 


3.0(0.1) 


0.3(0.1) 


4.5(0.3) 


G014.39-00.75A.. 


MMl 


17.84(0.04) 


1.0(0.1) 


1.2(0.1) 


0.8(0.4) 


0.4(0.1) 


2.9(0.4) 


0.1(0.1) 


3.0(0.9) 




MM2 


17.49(0.03) 


1.1(0.1) 


1.2(0.1) 


0.5(0.3) 


0.6(0.1) 


2.2(0.3) 


0.2(0.1) 


2.9(0.5) 


G014.39-00.75B.. 


MM3 


21.29(0.03) 


1.5(0.2) 


0.9(0.1) 


2.3(0.5) 


0.4(0.1) 


1.3(0.4) 






G014.63-00.57.. 


MMl 


18.75(0.01) 


4.0(0.4) 


1.8(0.1) 


2.2(0.1) 


2.4(0.1) 


2.4(0.1) 


1.2(0.1) 


2.7(0.1) 




MM2 


18.45(0.02) 


3.0(0.3) 


1.3(0.1) 


1.8(0.2) 


1.3(0.1) 


1.8(0.1) 


0.4(0.1) 


2.4(0.4) 




MM3 


17.64(0.04) 


1.0(0.2) 


1.4(0.1) 


2.1(0.4) 


0.5(0.1) 


1.9(0.2) 


0.2(0.1) 


1.5(0.3) 




MM4 


19.13(0.05) 


0.6(0.2) 


0.8(0.1) 


2.0(0.9) 


0.4(0.1) 


1.0(0.3) 






G016.93+00.24.. 


MMl 


23.80(0.01) 


1.5(0.1) 


0.9(0.1) 


1.7(0.2) 


0.5(0.1) 


1.7(0.3) 






G017.19+00.81.. 


MMl 


25.04(0.01) 


3.2(0.1) 


1.2(0.1) 


1.5(0.1) 


1.5(0.1) 


1.5(0.1) 


0.2(0.1) 


2.5(0.7) 




MM2 


22.75(0.01) 


3.1(0.1) 


1.3(0.1) 


1.3(0.1) 


1.7(0.2) 


1.6(0.1) 


0.5(0.1) 


2.3(0.3) 




MM3 


22.73(0.01) 


3.1(0.2) 


1.3(0.1) 


1.6(0.1) 


2.0(0.2) 


1.5(0.1) 


0.5(0.1) 


2.4(0.2) 




MM4 


21.64(0.06) 


0.5(0.1) 


2.8(0.1) 


0.9(0.3) 


0.3(0.1) 


3.6(0.3) 


0.2(0.1) 


2.6(0.4) 


G018.26-00.24.. 


MMl 


68.07(0.01) 


2.6(0.1) 


1.5(0.1) 


2.4(0.2) 


1.6(0.1) 


2.3(0.1) 


0.6(0.1) 


2.8(0.3) 




MM2 


67.75(0.02) 


2.6(0.1) 


2.0(0.1) 


2.8(0.1) 


1.7(0.2) 


2.6(0.2) 


0.3(0.1) 


2.8(0.3) 




MM3 


68.32(0.01) 


3.6(0.7) 


2.2(0.1) 


2.6(0.1) 


1.8(0.1) 


2.7(0.1) 


0.6(0.1) 


3.3(0.1) 




MM4 


68.31(0.02) 


3.4(0.4) 


2.1(0.1) 


2.2(0.1) 


1.8(0.1) 


2.8(0.1) 


0.8(0.1) 


2.9(0.2) 




MM5 


66.28(0.03) 


1.9(0.2) 


2.0(0.1) 


2.6(0.2) 


1.1(0.1) 


2.0(0.2) 


0.2(0.1) 


5.0(1.1) 


G022.06+00.21.. 


MMl 


51.17(0.02) 


1.8(0.1) 


1.7(0.1) 


1.9(0.2) 


1.5(0.1) 


2.2(0.2) 


0.5(0.1) 


3 6(^0 31 




MM2 


51.45(0.02) 


1.1(0.1) 


1.2(0.1) 


1.9(0.3) 


0.5(0.1) 


2.9(0.5) 


0.1(0.1) 


5.3(1.0) 


G023.38-00.12.. 


MMl 


98.42(0.02) 


1.7(0.1) 


1.6(0.1) 


2.3(0.2) 




1 Qcn 11 




3 MV\ 31 




MM2 


98.86(0.02) 


1.0(0.1) 


1.1(0.1) 


1.4(0.3) 


0.5(0.1) 


1.1(0.1) 






G024.37-00.15.. 


MMl 


58.83(0.03) 


0.9(0.1) 


2.1(0.1) 


2.8(0.3) 




9 (SiT) 91 


fl 9('fl 1 1 


A T(0 4.1 




MM2 


55.96(0.03) 


1.3(0.1) 


1.5(0.1) 


2.2(0.3) 


0.6(0.1) 


2.0(0.3) 


0.1(0.1) 


3.2(1.0) 


G024.61-00.33.. 


MMl 


42.80(0.02) 


1.3(0.1) 


1.4(0.1) 


1.2(0.2) 


fl f\<C\ 1 ~1 


1.5(0.2) 


fl 1 ("fl 1 1 


AQH 11 




MM2 


43.58(0.01) 


1.9(0.1) 


0.8(0.1) 


1.7(0.2) 


0.8(0.1) 


1.3(0.1) 


0.2(0.1) 


3 81 


G024.94-00.15.. 


MMl 


47.10(0.01) 


2.8(0.1) 


1.5(0.1) 


2.3(0.1) 


1.3(0.2) 


1.9(0.2) 


0.4(0.1) 


2.4(0.3) 




MM2 


48.03(0.02) 


1.8(0.1) 


1.4(0.1) 


2.2(0.2) 


0.8(0.1) 


2.1(0.2) 


0.2(0.1) 


3.7(0.7) 


G025.79+00.81.. 


MMl 


49.70(0.01) 


2.0(0.1) 


1.1(0.1) 


2.1(0.1) 


fl 6Cfl 1 1 


1 7^n 0\ 








MM2 


49.81(0.02) 


0.9(0.2) 


0.9(0.1) 


1.3(0.4) 










G030.90+00.00A.. 


MMl 


74.82(0.03) 


1.6(0.1) 


1.5(0.1) 


1.7(0.3) 


fl Qffl 1 1 


9 nm 1 1 


fl 3^ 1 1 


9 1 ("fl 4.1 


G030.90+00.00B.. 


MM2 


109.55(0.06) 


1.4(0.1) 


2.5(0.1) 


1.9(0.3) 










G030.90+00.00C.. 


MM3 


93.78(0.03) 


1.2(0.1) 


1.1(0.1) 


2.7(0.4) 


0.7(0.2) 


2.3(0.4) 


0.1(0.1) 


1 1(0 91 


G030.90+00.00D.. 


MM4 


37.61(0.05) 


0.8(0.1) 


1.0(0.1) 


1.9(0.8) 










G034.71-00.63.. 


MMl 


44.61(0.02) 


1.7(0.1) 


1.7(0.1) 


1.1(0.1) 


fl Sr'fl n 


2.2(0.2) 


fl 3Cfl 1 1 


4 4Cfl 31 




MM2 


45.43(0.02) 


1.9(0.4) 


1.3(0.1) 


3.0(0.3) 


0.6(0.1) 


1.5(0.1) 


0.2(0.1) 


6.4(0.8) 




MM3 


46.16(0.02) 


1.1(0.1) 


1.4(0.1) 


0.6(0.2) 


0.4(0.1) 


1.3(0.2) 






G034.77-00.81.. 


MMl 


43.36(0.08) 


0.6(0.1) 


1.6(0.2) 


0.5(0.7) 










G034.85+00.43.. 


MMl 


55.52(0.03) 


1.1(0.3) 


0.8(0.1) 


1.8(0.5) 


0.3(0.1) 


0.8(0.1) 








iVliVl 1 


jj.^jyyj.yjo) 


U. /{\J. 1 ) 


1 QfCi 1 ^ 


u. i\yj.^) 


0.3(0.1) 


2.2(0.4) 


0.2(0.1) 


4.1(0.6) 


G035.49-00.30B.. 


MM2 


45.31(0.01) 


2.7(0.3) 


1.0(0.1) 


3.5(0.3) 


0.8(0.1) 


1.3(0.1) 


0.1(0.1) 


2.4(1.2) 




MM3 


45.71(0.02) 


2.1(0.1) 


1.6(0.1) 


2.5(0.2) 


0.8(0.1) 


1.6(0.1) 


0.1(0.1) 


1.9(0.7) 


G037.44+00.14A.. 


MMl 


18.37(0.02) 


1.3(0.2) 


0.8(0.1) 


0.6(0.3) 


0.4(2.0) 


0.5(2.2) 






G037.44+00.14B.. 


MM2 


40.17(0.11) 


0.4(0.1) 


1.3(0.3) 


0.5(1.1) 










G050.06+00.06.. 


MMl 


54.11(0.03) 


0.9(0.1) 


1.3(0.1) 


1.5(0.3) 


0.3(0.1) 


1.3(0.3) 








MM2 


54.53(0.03) 


0.8(0.1) 


1.3(0.1) 


1.2(0.3) 


0.2(0.1) 


3.2(0.7) 






G053.81-00.00.. 


MMl 


24.13(0.02) 


1.3(0.1) 


1.4(0.1) 


1.3(0.2) 


0.3(0.1) 


3.1(0.9) 







